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FORJEWORD 


This  report  was  prepared  by  personnel  of  the  School  of  Engineering  of 
the  University  of  Dayton,  under  USAF  Contract  F33615-73-C-3023,  Project 
1369,  Task  01,  Work  Unit  17.  This  work  was  accomplished  under  the 
direction  of  the  Vehicle  Equipment  Division,  Air  Force  Flight  Dynamics 
Laboratory,  Air  Force  Wright  Aeronautical  Laboratories,  Wright-Patterson 
Air  Force  Base,  Ohio,  Mr.  George  J.  Sperry,  AFFDL/FEM,  Project 
Engineer.  ASD/PPMNA  was  the  controlling  procvxrement  office. 

This  report  covers  work  conducted  from  1 January  1974  to  15  May  1975. 
A prior  report,  AFFDL-TR-74-6,  entitled  "Landing  Gear/Soil  Interaction 
Development  of  Criteria  for  Aircraft  Operation  on  Soil  During  Turning  and 
High  Speed  Straight  Roll,  " was  also  developed  under  this  contract  and 
released  January  1974.  The  first  report  discusses  a turned  tire  test  program, 
initial  work  on  a digital  program  for  turning  prediction,  high  speed  drag 
prediction  and  a takeoff  program  development.  The  second  year's  effort, 
reported  herein,  discusses  the  use  of  the  turned  tire  data,  use  of  multipass 
data,  roughness  analysis  through  a digital  program,  and  development  of 
improved  predictive  digital  programs  for  takeoff,  landing,  turning  and  taxi 
performance.  The  authors'  submittal  date  was  1 October  1975. 

The  authors  wish  to  thank  Mr.  Sperry  for  his  support  and  assistance 
in  the  accomplishment  of  the  research  activities  and  in  the  successful 
conduction  of  these  training  sessions. 
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SECTION  1 

INTRODUCTION  AND  SUMMARY 

The  continuing  design  and  operational  requirement  for  aircraft  to 
operate  on  unimproved  runways,  as  evidenced  by  the  current  design 
requirements  on  the  MSTOL  aircraft  currently  under  development,  has 
resulted  in  continuing  efforts  to  identify  and  analyze  the  primary  and 
secondary  variables  which  influence  aircraft  flotation /ope ration  perform- 
ance on  soil  runways.  These  flotation /ope ration  variables  have  been  defined 
(1-51 

previously  . These  variables  include  drag,  sinkage,  braking,  turning, 
multipass,  high  speed,  soil  strength,  etc. 

The  current  research  effort  is  a part  of  a long  i-ange  research  program 
sponsored  by  the  United  States  Air  Force  Flight  Dynamics  Laboratory.  The 
objective  of  this  continuing  research  program  is:  (1)  to  identify  and  analyti- 
cally define  landing  gear/soil  interaction;  (2)  to  develop  criteria  for  establishing 
working  range  conditions  for  aircraft  in  their  landing,  takeoff,  braking,  and 
turning  modes  of  operation;  (3)  to  develop  systematic  design  procedures  for 

optimizing  the  flotation  and  surface  operations  capability  of  existing  and 

(5) 

future  aircraft.  Phase  I of  this  program  included  a survey  of  the  flotation 

problem,  establishment  of  the  critical  parameters,  and  an  investigation  of 

available  flotation  data  leading  to  the  development  of  a flotation  analysis 
(4) 

equation.  Phase  II  included  the  development  of  an  empirical  sinkage 

prediction  equation,  development  of  a lumped  parameter  simulation  prediction 

technique,  conducting  the  Rolling  Single  Wheel  Verification  Tests,  and  the 

development  of  the  Single  Wheel  Relative  Merit  Index  (RMI)  system  for 

(3) 

defining  comparative  flotation  capacity.  Phase  III  - Part  I consisted  of 
the  development  of  the  multiwheel  sinkage -drag  analysis  equations,  conducting 
the  Multiwheel  Verification  Tests,  and  the  development  of  a lumped  parameter 
iteration  technique  for  simulating  the  interaction  of  dual  tires  on  soil. 


Phase  III  - Part  II  ^ ^ included  the  Braked  Wheel  Verification  Tests,  the 
development  of  a lumped  parameter  braking  simulation  technique  computer 
program,  the  development  of  the  braking  analysis  equations  for  defining 
braking  drag  ratios,  and  preliminary  studies  of  multipass  and  high  speed 
effects  of  aircraft  tire  operation  on  soil. 

The  current  two-year  effort  was  aimed  at  developing  analysis 
techniques,  conducting  verification  tests,  and  developing  predictive 
equations  for  landing  gear/soil  interaction  as  necessary  to  estimate 
vehicle  performance  during  turning,  high  speed  and  multipass  operations. 

During  the  first  year  the  following  tasks  were  completed. 

1)  A turned  tire  test  program  was  conducted  at  UD  and  AEWES 
to  provide  measured  data. 

2)  An  Aircraft  Turning  Operation  Predictive  Computer  Program 
was  evolved. 

3J  High  speed  drag  ratio  analysis  and  predictive  techniques  were 
generated. 

4)  An  existing  takeoff  length  program  was  modified. 

5)  A Landing  Gear/Soii  Interaction  Training  Session  was 
conducted. 

During  the  second  year  the  following  have  been  accomplished. 

1)  The  Aircraft  Turning  Program  was  completed  and  had  incorporated 
into  it  the  complete  results  of  the  turned  tire  test  data.  Additionally, 
the  routine  now  has  the  capability  to  accept  thrust  differential  as 
well  as  differential  braking  and  nose  wheel  steering.  Examples 

of  the  response  of  three  vehicles  are  presented. 

2)  Multipass  data  were  collected  at  UD  and  AEWES  to  provide 
measured  data  for  single  wheels,  rolling  and  braked,  over 
alternating  paths. 
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3)  Predictive  equations  were  evolved  for  single  wheel  multipass 
operations. 

4)  A landing  length  program  was  developed  by  modifying  an  existing 
routine. 

5)  Existing  tire/soil/roughness  data  were  reviewed  and  an  analytical 
model  developed  to  predict  soft  tire/soft  soil  response.  The 
digital  program  was  completed  but  not  verified  using  test  data. 

6)  Existing  start-up  force  data  were  examined  to  establish  preliminary 
criteria  for  estimating  start-up  forces  in  relation  to  free  rolling 
drag. 

7)  A master  program  was  developed  to  permit  a user  to  calculate 
takeoff  distance,  landing  distance,  turning  response,  and  multi- 
pass performance  for  one  set  of  data. 

8)  A test  plan  was  developed  to  indicate  how  testing  could  be  conducted 
to  verify  the  computer  generated  data  of  the  developed  digital 
programs. 

9)  A second  Landing  Gear/Soil  Interaction  Training  Session  was 
conducted. 

The  research  conducted  has  provided  additional  capability  in  determining 
tire/soil  performance  for  all  aspects  of  unprepared  runway  operations.  For  a 
particular  mission,  the  vehicle  lands,  rolls  out,  turns,  taxiis,  on^loads  or 
off-loads,  taxiis  back  out,  turns,  and  then  takes  off  perhaps  to  repeat  the 
cycle.  The  results  of  this  program  and  others  provide  data  to 

estimate  the  landing  gear  performance  for  these  mission  segments  in  terms 
of  drag,  sinkage  and  rut  depth. 
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SECTION  II 

TURNED  TIRE /SOIL  INTERACTION  ANALYSIS  AND  APPLICATION 


A.  INTRODUCTION 

Aircraft  turning  operation  is  an  import^.nt  consideration  in  the  flotation 
analysis  of  aircraft  operating  on  unprepared  runways.  Landing  gear  failures 
due  to  high  side  loads  brought  about  by  turning  either  in  a high  sinkage  con- 
dition or  out  of  accumulated  ruts  are  possible.  The  forces  developed  can 
create  large  steering  torque  requirements.  Additionally,  the  turning  process, 
whether  by  steering  or  braking,  contributes  to  the  increased  deterioration 
of  the  runway  surface  making  further  operations  questionable. 

The  uncertainty  associated  with  a turning  operation  is  reflected  in 
current  design  and  analysis  criteria.  The  C-130  operations  manual  indicates 
that  there  are  restrictions  on  the  steering  angle  as  a function  of  taxi  speed, 
but  not  soil  type,  or  steering  rate.  Military  Specification  MIL-A-008862A, 
Airplane  Strength  and  Rigidity,  indicates  <-hat  there  are  side  force,  drag 
force  and  lateral  force  ratios  applicable  to  operations  in  soils,  but  does  not 
mention  the  effects  of  steering  angle,  steering  rate,  velocity,  soil  type  or 
braking  upon  turning. 

There  are  several  means  available  to  turn  an  aircraft.  The  nose  wheel 
cam  be  turned  at  an  angle  relative  to  the  aircraft  plane  of  symmetry  and  the 
side  forces  developed  at  the  nose  tire  will  turn  the  aircraft.  The  engines  can 
be  used  to  turn  the  vehicle  by  using  differential  thrust.  Asymmetrical  braking 
creates  a moment  that  turns  the  vehicle.  Any  combination  of  the  nose  wheel 
steering,  braking  and  thrust  differential  cam  also  be  used.  Regardless  of 
the  meams  used,  all  tires  will  have  some  turn  angle  which  generates  a side 
force.  This  is  shown  in  Figure  1 which  represents  the  geometry  of  a vehicle 
turning  due  to  some  applied  moment.  Notice  that  all  wheels  have  some  turn 
angle.  Frequently  in  calculations  of  turn  radii  figures  will  be  drawn  showing 
the  turning  aixis  as  the  intercept  of  a line  through  the  main  gear  and  one 
perpendicular  to  the  plane  of  the  nose  wheel.  This  is  obviously  erroneous 
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Figure  1.  Aircraft  Turning  Geometry, 
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for  nose  wheel  steering  only  turns, since  if  all  three  wheels  have  zero  turn 
angle,  there  is  no  applied  moment  to  turn  the  vehicle.  As  will  be  discussed 
later,  the  true  turning  of  the  vehicle  requires  an  equilibrium  condition  of 
side  and  drag  forces  at  the  wheels,  all  at  different  turn  angles,  with  the 
centrifugal  forces  in  the  turn.  Simultaneously  the  moments  generated  by 
tliese  factors,  brakes  and  engines  must  be  properly  accounted  for. 

For  computation  of  the  turning  response  it  is  necessary  to  know  how 
the  side  and  drag  forces  of  the  tire  are  influenced  by  the  turn  angle.  A 
review  of  the  existing  literature  indicated  very  little  was  known  analytically 
or  experimentally  about  side  loads  on  aircraft  tires  turning  on  pavement, 
with  less  known  about  turning  in  soils.  And  yet.  Figure  1 indicates  that 
each  wheel  will  have  a turn  angle  whether  dictated  by  the  turning  radius,  as 
for  the  main  gear,  or  as  by  the  oifference  between  steering  angle  and  turning 
angle  at  the  nose.  The  information  required  was  found  by  conducting  turned 
tire  tests  with  aircraft  tires  towed  at  a constant  turn  angle,  0,  down  a test 
tr-'-'k. 


B.  TURNED  TIRE  TEST  PROGRAMS 


Two  series  of  turned  tire/soil  interaction  tests  were  conducted.  The 
purpose  of  both  test  programs  was  to  obtain  experimental  data  to  permit 
further  studies  of  the  variables  that  influence  the  performance  of  aircraft 
tires  while  operating  in  a turned  (yawed)  mode.  The  test  program  using 
the  Universit/  of  Dayton  Linear  Tire/Soil  Test  Track  (see  Figure  2 ) was 
conducted  first  (March  1973)  and  the  results  were  used  primarily  to  design 
the  test  program  conducted  at  the  Army  Engineers  Waterways  Experiment 
Station  (AEWES)  (see  Figures)  small  scale  facility  during  June  and  July 
1973,  The  specific  objectives  of  the  turned  tire  tests  included: 

1.  Establishment  of  relationships  between  in-line  drag  (see  Figure  4) 
side  load,  and  lateral  load  as  a function  of  turn  angle. 


rogram  using  University  of  Dayton  Linear  Tire/Soil  Test  Track 
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2.  Determination  of  the  variations  in  percent  slip  (negative)  with 
increasing  turn  angle. 


3.  Establishment  of  the  influence  of  increasing  turn  angle  on  the 
resultant  sinkage  and  rut  depth. 


The  specific  test  results  discussed  below  are  contained  in  detail  in 


References  13  and  14. 


Soil  Type,  Test  Equipment  and  Sinkage  Ranges 


The  AEWES  tests  were  conducted  in  mortar  sand  and  buckshot  clay. 
The  soil  strength  varied  from  cone  index  readings  of  approximately  35  to 
nearly  45  in  sand,  and  37  to  74  in  clay.  Applied  vertical  loads  varied  from 
550  to  1,  030  pounds  for  the  sand,  and  1,  310  to  1,  580  pounds  for  the  clay. 
These  values,  along  with  35  to  40  percent  deflections  on  8:50-10  tires  and 
7:00-6  tires;  provide  a 0°  turn  sinkage  ratio  of  0.03  to  0.04.  The  sinkage 
ratio  is  defined  as  the  ratio  of  sinkage  Z,  to  tire  diameter  D,  and  the  values 
selected  are  indicative  of  operational  conditions  on  soil  runways. 


The  University  of  Dayton  tests  were  conducted  on  riverwash  sand  and 
clayey  silt.  Soil  strengths  of  30  to  46  cone  index  values  for  clayey  silt  and 
70  to  80  for  riverwash  sand  were  selected.  These,  with  applied  loads  of 
115  to  145  povmds  on  a 5:00-5  Type  III  tire  at  17  percent  deflection,  generate 
the  desired  sinkage  ratio. 


Test  Parameters  Measured 


For  both  series  the  following  were  measured. 


vertical  load,  P 


drag  load  (in-line  drag),  R 

side  load  (perpendicular  to  plane  of  wheel),  S 

resultant  lateral  force  (determined  from  R and  S above),  L 


horizontal  carriage  velocity 


turn  angle,  0 


tire  deflection,  d 

tire  contact  geometry  (rigid  surface  contact  area),  A 

soil  properties  (moisture,  density,  classification  properties. 

Cone  index.  Cl  ) 
avg 

permanent  rut  depth  (instantaneous  sinkages  were  determined  from 
rut  depth  and  rebound  measurements) 

All  University  of  Dayton  tests  were  conducted  at  six  feet  per  second. 

All  AEWES  data  were  collected  at  ten  feet  per  second. 

C.  TEST  RESULTS 
Riverwash  Sand 

The  curves  shown  in  Figure  5 summarize  a total  of  22  tests 
conducted  on  riverwash  sand.  The  buildup  of  side  load  and  resultant  lateral 
load  with  increased  turn  angle  is  evident. 

Clayey  Silt 

The  results  of  23  tests  in  clayey  silt  are  shown  in  Figure  6 and  indicate 
the  significant  increase  in  side  force  with  increasing  turn  angle.  The  in-line 
drag  is  relatively  constant  over  the  turn  angles  measured. 

Mortar  Sand 

Twenty-one  tests  in  mortar  sand  generated  curves  very  similar  to 
those  of  riverwash  sand  as  shown  in  Figure  7 . Very  little  difference  exists 
between  the  two  sand  test  series  although  two  different  size  tires  were  used. 

Buckshot  Clay 

The  results  of  23  tests  in  buckshot  clay  are  shown  in  Figure  8 • The 
tests  conducted  with  larger  tires  in  buckshot  clay  indicate  less  side  force 
for  a given  turn  angle  than  that  measured  in  riverwash  sand  using  a smaller 
tire. 


11 


10*  15* 

Turn  Angle,  d> 


Figure  8.  AEWES  Test  Resvilts  for  Clay. 


Comparative  Summary  Analysis  ot'  Turned  Tire  Tests 


The  results  of  both  turned  tire  test  programs  were  analyzed  and 
compared.  The  results  of  this  comparison  indicated  the  following  points 
of  similarity. 

1.  Sinkage  and  rut  depth  increase  markedly  with  increasing  turn 
angle  for  aircraft  tires  operating  on  sand  (frictional)  type  soils, 
but  only  moderate  increases  in  sinkage  occur  on  clay  type  soils. 
See  Figures  9 and  10  . 

o 

2.  The  in-line  drag  ratio  remains  relatively  constant  over  the  0 
to  20°  turn  angle  range  for  aircraft  tires  on  clay  tyre  soils. 

3.  Significant  amounts  of  soil  buildup,  opposite  to  the  direction  of 
turn,  occur  for  turning  tires  in  both  sand  and  clay  types  of  soil. 

Perhaps  the  most  significant  result  of  the  turned  tire  tests  is  some 
preliminary  indication  of  the  magnitudes  of  side  load  buildup  on  a turned 
tire  operating  on  a soil  runway.  Figure  11  presents  this  type  of  result  as 
related  to  the  nonturned  tire  sinkage  ratio  (Z,  x3).  Reference  to  Figure  11 
shows  that  the  magnitude  of  this  sinkage  has  a considerable  influence  on  the 
resulting  lateral  loads.  Lateral  load  ratios  may  very  well  exceed  one-half 
in  situations  of  tight  turn  angles  and  high  sinkage  conditions.  Some  of  the 
general  conclusions  as  to  turned  tire  performance  on  soil  runways  would 
include; 

1.  The  lateral  load  ratio  (L/P)  very  likely  will  approach  and  in  some 
cases  exceed  one -half  for  aircraft  operating  on  moderate  strength 
soil  runways  (rolling  drag  ratios  greater  than  0.  10)  where  tight 
turning  situations  are  encountered. 

2.  The  tire  negative  slip  will  increase  throughout  the  entire  turn 
angle  range  of  0°  to  20°.  Negative  slips  approached  25%  in  sand 
and  15%  for  clay  in  the  20°  turn  situation. 
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Figure  9.  Turned  Tire  Rut  Profile,  Clay,  (AEWES). 
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Figure  10.  Turned  Tire  Rut  Profile,  Sand,(  AEWES). 
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3.  Unlike  the  behavior  of  turned  tires  on  rigid  pavement,  the  percent 
tire  deflection  had  little  influence  on  the  lateral  load  ratio  for 
either  sand  or  clay  operation. 

4.  More  severe  rutting  will  occur  for  aircraft  turning  operations 
on  frictional  type  soils  (sand)  than  on  cohesive  (clay)  type  soils. 

5.  No  clear  trend  existed  in  the  variation  of  pneumatic  trail  (r)  with 
turn  angle.  The  pneumatic  trail  increased  with  increasing  turn 
angle  in  sand  but  decreased  with  increasing  turn  angle  in  clay. 

Although  sufficient  turned  tire  tests  were  not  conducted  to  encompass 
a broad  range  of  sinkage  ratios  (Z/D),  the  results  do  indicate  that  an 
increase  in  the  sinkage  ratio  v/ill  result  in  an  incieased  lateral  load  ratio 
for  a turned  tire  for  all  other  conditions  being  constant. 


D.  AN  AIRCRAFT  TURNING  MANEUVER 


The  side  force  data  collected  were  incorporated  within  a digital  program 
to  compute  the  turning  performance  of  a tricycle  landing  gear  aircraft.  The 
physical  description  of  the  vehicle  in  a turn  is  shown  in  Figure  12  . The 
symbols  are  defined  as: 


0 

0 

0 


n 

mjt 

mr 


nose  wheel  turning  angle 

left  main  gear  wheel  turning  cingle 

right  main  gear  wheel  turning  angle 


W 

R 


V 


weight  of  aircraft 

drag  loads  (braking  or  driving)  in  longitudinal  wheel 


direction  R^.  R^.  R^^) 


side  load  perpendicular  to  wheel  longitudinal  direction 


(S  ,,  S ,,  S S and  S .) 
nl  ml  m2  m3  m4 


aircraft  forward  velocity 
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Q = instantaneous  radius  of  turn 

B = caster  angle  of  nose  wheel 

E = distance  between  main  mounts 

AVM  = main  gear  center  of  gravity  to  average  vehicle  center  of 

gravity 

AL  = nose  gear  center  of  gravity  to  average  vehicle  center  of 
gravity 

g = gravitational  constant 


The  aircraft  at  any  instant  is  turning  about  some  instantaneous  center 
of  rotation.  The  side  forces  are  generated  by  wheels  being  forced  to  follow 
a circular  arc  not  in  line  with  the  plane  of  the  tire.  The  ma\n  wheels  do  so 
because  of  the  radius  of  curvature  generated  by  all  forces,  whereas  the 
nose  wheel  side  force  is  generated  by  both  the  radius  of  turn  and  the  steering 
angle.  As  the  vehicle  travels  about  the  turn  at  a given  velocivy,  the  sunnmation 
of  side  and  drag  forces  create  the  net  force  which  generates  a centripetal 
acceleration.  The  forces  generate  moments  about  the  center  of  gravity  of  the 
vehicle  which  must  be  proportional  to  the  angular  acceleration. 

A summation  of  forces  and  moment  about  the  center  of  gravity  yields; 

S Fy  = M’y  = LS  + S cosB  + R sinB 

m n n 


^ ^CG  = ' e = ' 


^ V4>  > f 


(S  ^+S  J)xAVM+S  (cosB)AL-R  (sinB)  AL. 
m3  m4  n u 

Over  a small  increment  of  time  if  constant  velocity  motion  is  about  a fixed 
radius,  the  radial  acceleration  of  the  body  is  that  of  centripetal  acceleration, 
and  the  angular  acceleration  is  dictated  by  the  change  in  angular  velocity 
from  one  radius  to  the  next. 

vSince  the  centrifugal  force  is  known  for  any  assumed  radius. 


I 


19 


centrifugal  force  = -r — . The  inertial  force  and  moments  can  be  placed 

R 

on  the  right  side  of  the  equation  and  then  both  the  equations  are  of  the  form, 

LF  = 0 

S M = 0 . 

The  problem  is  to  find  forces  which  will  simultaneously  satisfy  both  equations. 

A Newton -Raphson  technique  is  used  to  iterate  from  an  assumed 
radius  of  curvature  to  one  that  provides  equilibrium  of  applied  and  inertial 
forces  and  moments.  Other  subroutines  calculate  the  weight  distribution, 
sinkage,  drag,  turn  angles,  multiwheel  effects,  and  thrust  required. 

The  program  accepts  nose  wheel  steering  angle  as  a function  of  time 
along  with  wheel  slip  as  a function  of  time.  These,  along  with  a complete 
description  of  the  vehicle,  landing  gear  and  tire  configurations,  are  required 
as  input.  The  program  then  calculates  the  path  over  the  ground,  the  tire/ 
soil  interaction  forces,  and  the  vehicle  acceleration. 

A complete  description  of  !±ie  routine  developed  during  the  first  year 
of  the  research  program  is  available  in  Reference  15.  All  subroutines 
were  fully  discussed  and  sample  results  presented.  At  the  time  the  digital 
program  had  several  short-comings  and  these  were  mentioned  in  the  con- 
clusions of  the  referenced  report.  Specifically, 

1.  The  braked  drag  equation  was  an  approximation  to  other  equations 
that  were  available. 

2.  Pneumatic  trail  data  were  not  included. 

3.  Throttle  controls  were  not  possible. 

4.  Corrections  were  not  included  for  "driven"  and  "pulled"  wheel. 

5.  Direction  reversals  were  not  possible. 
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Additionally,  the  program  was  restricted  to  CBR  6 clay  soil.  Sample 
computations  were  made  using  the  initial  turning  program.  The  C-130E  was 
chosen  as  a vehicle  that  could  Se  analyzed  by  the  developed  program.  Most 
of  the  aircraft  input  data  required  were  available  and  some  information  from 
the  pilots  handbook  provided  an  insight  into  how  the  aircraft  is  turned. 
Specifically,  the  pilots  handbook  states  that  60°  nose  wheel  steering  is 
permitted  up  to  5 knots.  Nose  wheel  steering  of  20°  is  permitted  at  taxi 
speeds  of  up  to  20  knots.  For  these  conditions,  what  are  the  turn  radii,  side 
forces  and  drag  forces,  and  how  do  they  compare  with  the  criteria  of 
existing  military  standards? 

The  first  case,  60°  nose  wheel  steering  at  5 knots,  was  examined  by 
assuming  a time-steering  angle  relation.  It  was  assumed  that  it  would 
take  10  seconds  to  reach  a 60°  steering  angle.  The  vehicle  will  turn  at  an 
equilibrium  radius  of  51  feet  and  require  a minimum  runway  width  of  106 
feet.  The  width  includes  the  entry  transient  as  well  as  the  landing  gear 
width.  The  centrifugal  force  generates  an  aircraft  side  load  factor  of  0.054, 
and  thrust  required  reaches  a maximum  of  23,000  pounds,  which  is  a signifi- 
cant portion  of  the  40,000  pounds  available.  The  response  is  shown  in 
Figure  13. 

The  force  ratios  developed  during  the  turn  indicate  that  while  the 
main  gear  develops  very  little  side  force,  the  nose  wheel  reaches  the  maxi- 
mum pernriissible.  Therefore,  the  60°  nose  wheel  steering  maneuver  would 
appear  to  be  that  which  is  limited  by  the  side  force  capability  of  the  nose 
gear. 

The  second  case,  that  of  20°  nose  wheel  steering  is  examined  by 
beginning  with  a 5 knot  taxi  speed  and  then  increasing  it  to  10  and  20  knots. 

For  these  turns  it  was  assumed  that  3-1/2  seconds  v/ere  required  to  turn 

the  nose  wheel  to  20°.  The  computed  results  indicate  that  in  going  from  5 

knots  to  20  knots,  the  minimum  runway  width  goes  from  214  to  394  feet, 

and  the  aircraft  side  load  factor  from  0,22  to  0,40.  The  force  ratios  indicate  that 

again  the  main  gear  does  not  develop  the  allowable  levels  but  that  the  nose 

wheel  does. 
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Feet  of  Rxmway  Length 


Figure  13.  Nose  Wheel  Steering  at  60°, 


As  another  case  to  be  examined,  it  was  assumed  that  both  braking  and 
nose  wheel  steering  would  be  used.  Braking  at  32  percent  slip  on  the  inboard 
gear  and  sufficient  nose  wheel  steering  to  "trail"  the  turn  were  used.  The 
braked  turn  at  5 knots  generates  an  85  foot  minimum  runway  width  with  a 
vehicle  side  force  ratio  of  0.062.  The  landing  gear  side  force  ratios  are 
again  quite  s/nall  although  the  drag  forces  developed  were  limiting.  It  is 
interesting  to  note  that  the  turn  computed  is  barely  possible,  since  the 
thrust  required,  37,700  pounds,  approaches  that  available  for  the  C-130E 
aircraft.  See  Figure  14  . 

Improved  Turning  Program 

During  the  second  year  several  goals  were  established  to  improve  the 
existing  routine.  These  were: 

1)  Incorporate  sinkage  and  drag  relations  from  sand  and  clay 

soil, 

2)  Incorporate  new  braked  drag  relations, 

3)  Utilize  turned  tire  test  results  to  incorporate 

a)  pneumatic  trail  for  sand  and  clay 

b)  side  force  ratios  for  sand  and  clay 

c)  slip  for  sand  and  clay  as  a function  of  turn  angle 

d)  in-line  drag  modifications  for  sand  and  clay 

4)  Add  thrust  variations  with  time, 

5)  Add  the  longitudinal  degree  of  freedom. 

All  of  these  were  not  achieved,  but  significant  advancements  were 
made.  In  the  following  paragraphs  the  raw  data  from  which  the  equations 
were  evolved  may  not  be  presented.  The  complete  information  if  required 
is  contained  in  Reference  13. 
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Pneumatic  Trail 


Pneumatic  trail  (r)  data  were  plotted  for  the  8:50  and  7:00  tire  for  both 
clay  and  sand.  For  clay  the  trail  versus  turn  angle  was  plotted,  then  trail 
over  footprint  length  (1)  versus  turn  angle  to  be  consistent  with  published 
literature  (16).  The  wide  scatter  '.n  the  data  indicated  that  some  other  term 
would  be  necessary  to  combine  the  data  into  some  easily  usable  form.  Not 
surprisingly,  instantaneous  sinkage  (Z)  was  found  to  pull  the  data  together 
in  terms  of  tire  deflections  as  shown  in  Figure  15  for  clay  soil. 

From  the  figure  it  is  seen  that  pneumatic  trail  is  behind  the  aixis  of 

J* 

the  tire  by  the  negative  of  — . Z,  and  that  the  trail  decreases  with 
increased  turn  angle.  Secondly,  the  curves  are  grouped  according  to  tire 
deflection  with  the  40  percent  deflection  tires  having  greater  trail  values 
for  a given  turn  angle. 

Sa,-cral  approximate  expressions  for  the  trail  were  attempted.  l inally, 
the  simplest  expression  evolved  was 


r X 


Z 

I 


0.00675 
sin  8 


(2.0) 


where 


r is  tlie  trail  (inches) 

Z is  the  instantaneous  sinkage  (inches) 

I is  the  footprint  length  (inches) 

8 is  the  turn  angle 

DE  is  the  deflection  in  percent. 


The  equation  provides  a smooth  transition  from  30  to  40  percent 
deflection  which  was  assumed  to  be  sufficient  for  most  tires  of  interest.  The 
equations  specifically  "fits"  the  35  percent  and  40  percent  deflection  data 
measured. 


Similar  procedures  were  followed  for  sand.  The  data  aro  shown  in 
Figure  16.  The  curves  are  entirely  different  in  that  the  trail  is  positive, 
ahead  of  the  tire  axis,  and  increases  with  turn  angle.  The  curves  are 
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0 at  40%  Deflection 

8:50  at  40% 

r-  Predictive  Equati< 


Turn  Angle  (Degrees) 


Figure  16.  Pneumatic  Trail  Parameter  versus 
Turn  Angle  for  Cohesionless  Soil. 


functions  of  both  deflection  and  sinkage.  Equations  developed  were; 


r X 


= 0.0375.  10  p ^ Z s 0.6" 


where 


where 


P is  the  turn  angle  in  degrees,  and 

Z > 0.6", 

DE  . 


r X - = 0.  1075.  10"^  P K 

I 


K = 4.5 


10 


in  order  to  correct  for  the  influence  of 


deflection  at  greater  sinkages. 


The  equations  listed  above  were  incorporated  into  the  turning  program. 

During  the  incorporation  of  the  pneumatic  trail  into  the  program, 
calculations  were  made  for  a small  fighter  type  aircraft  of  approximately 
30,  000  potinds.  For  single  wheel  landing  gear  in  sand  and  clay  it  was  found 
that  the  pneumatic  trail  made  negligible  difference  in  the  path  over  the 
ground  or  thrust  required.  This  seems  reasonable  in  that  the  torque 
generated  at  the  tire  is  very  small  in  proportion  to  the  moment  generated 
by  the  side  force  about  the  aircrt^t  center  of  gravity.  Therefore,  the 
torque  generated  by  the  inclusion  of  pneumatic  trail  is  important  from  the 
standpoint  of  steering  control  torque  values,  but  unimportant  in  trajectory 
calculations.  Since  the  torque  for  one  tire  is  always  the  side  force  times 
the  trail  and  the  rigid  body  moment  is  always  the  side  force  times  the 
distance  to  the  center  of  gravity  of  the  vehicle;  it  seems  reasonable  to 
assume  the  relative  importance  of  the  trail  established  would  be  the  sam.e 
for  multiple  wheel  configurations. 


Side  Force  Ratios 

Side  force  ratio  data  for  sand  euid  clay  are  shown  plotted  in  Figures 
17  and  18.  The  curves  show  little  dependence  upon  tire  deflection  or 
size  and  hence  two  equations  were  found  to  fit  the  data.  These  were; 
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e Force  Ratio  (S/P) 


s 

p 


0.438  sin  (4p)  for  clay  and 


P 


0.  432  sin  (48)  for  sand 


Slip  Data 

Free  rolling  slip  generated  by  the  turned  tire  is  plotted  in  Figures 
19  and  20.  No  attempt  was  made  to  establish  sinkage  or  tire  deflection  of 
tire  size  dependence.  Approximate  curves  evolved  were: 

2 

S = 5. 00  + 0.0425  8 for  sand  and 

S = 3.  00  + 0,0325  8 for  clay 


where  8 is  in  degrees. 


The  slip  measured  was  used  in  the  braked  drag  equation  to  determine 
whether  or  not  drag  as  only  a function  of  slip  co\ild  be  computed  properly. 
This  did  not  work.  Drag  created  by  braking  in  terms  of  slip  can  be  cal- 
culated from  a braked  drag  relation.  Drag  increases  with  turn  angle  are 
apparently  not  functions  of  the  slip  due  to  the  turn  angle.  It  was  found, 
however,  that  calculated  instantaneous  sinkage  for  a turn  angle  can  be 
computed  using  the  relation 


Z 


3 

TsT 


which  is  similar  to  a portion  of  the  braked  drag  equation.  In  this  equation 

the  turned  instantaneous  sinkage,  Z , is  calculated  using  tlie  slip,  S,  due  to 

P 

the  turn  angle  and  the  maximum  fully  braked  sinkage  The  instantaneous 

sinkage  is  then  used  in  Ceilculating  rut  depth. 


In-Line  Drag 

As  mentioned  above,  the  in-line  drag  does  not  change  with  slip  as 
originally  hoped.  Consequently,  the  in-line  drag  ratios  were  plotted  as  shown 
in  Figures  21  and  22.  For  each  curve  it  was  realized  that  we  need  not  meet 
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In  Line  Drag  Ratio  (R/P)  ^ Line  Drag  Ratio  (R/P) 


8:50  at  35% 
8:50  at  40% 


IN  LINE  DRAG 
RATIO  VERSUS 
TURN  ANGLE 
IN  COHESION  LESS 
SOIL 


(AEWES) 


7:00  at  40% 
-7:00  at  35% 


Predictive  Equation 


Turn  Angle 


Figure  21.  In-Line  Drag  Ratio  versus  Turn  Angle  in 
Cohesionless  Soil. 


Predictive  Equation- 


8:50  at  40% 
■8:50  at  35% 


IN  LINE  DRAG  RATIO 
VERSUS  TURN  ANGLE 
IN  COHESIVE  SOIL 


(AEWES) 


7:00  at  35% 
7:00  at  40% 


Turn  Angle 


Figure  22.  In-Line  Drag  Ratio  versus  Turn  Angle  In 
Cohesive  Soil. 


the  magnitudes  shown,  but  relations  that  establish  changes  .‘’’•om  a zero  turn 
condition.  These  relations  were  found  and  were; 


( -) 

' P'  t 


R 


-2^  2, 


( - ) (1  +0.27  .10  3 ) for  sand 

* r\ 


{1  +0.25.  lO’^P)  for  clay 


R 


where  ^ "p  \ turned  drag  ratio 


(“  )^  is  the  rolling  drag  ratio  at  zero  turn  auigle, 

and  P is  the  turn  angle  in  degrees. 

By  defining  the  turned  drag  ratio  in  this  manner  it  is  possrble  to  calcu- 
late the  free  rolT.ng  drag  ratio  as  previously  found  and  correct  for  the  given 
turn  angle. 

Thrust  Control 

Thrust  control  was  included  by  incorporating  a constant  thrust  diffe'ren 
tial  in  order  to  maintain  a specified  taxi  speed.  If  thrust  is  controlled  on 
both  engines  as  a function  of  time,  what  are  the  relations  used  as  input? 

This  type  of  question  was  similar  to  that  asked  for  steering  and  braking 
control.  Given  that  the  program  user  wants  to  vary  the  thrust  or  the  braJcing, 
what  time  relation  is  used?  In  order  to  eliminate  some  of  this  problem  and 
yet  provide  useful  data,  it  was  decided  that  a specified  differential  of  power 
would  be  used.  That  is,  if  a differential  of  1.  00  is  specified,  the  left  engine 
is  1.00  times  the  right  engine  for  a symmetrical  thrust  condition.  If  0,  50 
is  specified,  the  left  engine /engines  generate  one -half  that  of  the  outboard 
engine.  The  program  calculates  the  thrust  required  to  maintain  a given  taxi 
speed.  This  value  is  then  used  to  calculate  the  thrust  on  each  side  and  a 
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calculated  moment  applied  to  the  summation  of  moments  equation.  The 
generation  of  differential  thrust  in  this  manner  is  hoped  to  simulate  the 
response  cf  a pilot  who  when  entering  a turn  wants  thrust  differential  while 
maintaining  the  same  taxi  speed. 

Braking  and  Rolling  Drag  Relations 

The  braked  and  free  rolling  drag  and  sinkage  equations  of  Reference 
17  were  included  in  the  digital  program.  Therefore,  the  drag  subroutines 
were  updated  to  permit  the  drag  to  be  calculated  for  both  cohesive  and 
cohesionless  soils  defined  hi  terms  of  their  cone  index  values.  Braked  drag 
relations,  as  well  as  free  rolling  drag,  were  available  for  both  soil  types. 

Additional  Revisions  Possible 

Corrections  for  "driven"  and  "pulled"  wheels  on  the  same  axle  were 
not  incorporated.  The  possibility  of  a turning  reversal  is  not  available. 
Lastly,  rotating  effects  of  wheels  and  engines  were  not  included.  All  three 
aspects  can  be  incorporated  later. 

E.  SAMPLE  CALCULATIONS 

The  resulting  improved  digital  program  can  best  be  described  by 
presenting  some  sample  calculations.  The  C-130  was  again  used  as  a 
representative  vehicle.  The  input  required  is  presented  in  Figure  23  in 
which  a computer  generated  list  appears  as  the  heading  of  each  complete 
turn  calculation. 

Most  of  the  data  listing  are  self  explanatory.  All  dimensions  not 
specifically  mentioned  are  in  inches.  The  cone  index  value  shown,  192.  00, 
corresponds  to  a CBR  value  of  6 since  32  to  1 was  used  for  cohesive  soil. 
Cohesive  soil  is  specified  by  NTYP  equal  1.  The  engine  thrust  parameter  of 
1.00  indicates  that  the  left  engine  develops  1.00  times  the  right  engine  for 
symmetrical  thrust  distribution.  Maximum  steering  angle  for  the  nose 
wheel  is  specified  at  60.00  degrees  and  the  meocimum  slip  is  0.00  percent. 
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ATRCOAPT  PARAMETERS 


GPOSS  WEIGHT 

GW  = 

1330CC. 

C LBS 

> • 

INITIAL  ground 

VEL 

. VO  = 

5 

.C  KNOTS 

TIRE  DATA  — 

NOSE 

MAIN 

TIRE  DEFLECTIONr/.) 

OE  = 

35. CC 

DEM 

- 

35. CO 

TIRE  DIAMETER 

(IN. 

) DN  = 

36. OC 

DM 

53.20 

FLANGE  OTAM. 

(IN. 

) OFN  = 

16.50 

OFM 

s 

24. CD 

TOTAL  NO.  OF  TIRES 

NN  = 

2 

NM 

- 

4 

NO.  OF  TIRES/BOGIE 

N1 

2 

NO.  OF  TIRES  — 

TANDEM 

NMl 

2 

NO*  OF  IIRES— 

TWIN 

NNl 

0 

TIRE  TYPE 

TYPE  3 

TYPE 

3 

TIOE  AND  C.G.  OISTANCES--(IN  INCHES) 

NOSE  GEAR  C.G. 

TO 

FWO  C.G. 

L 

2 

346.33 

NOSE  GEAR  C.G. 

TO 

AVE  C.G. 

AL 

s 

358,80 

NOSE  GEAR  C.G. 

TO 

AFT  C.G. 

LL 

r 

371.30 

NOSE  GEAR  C.G. 

TO 

MAIN  GEAR 

C.G. 

F 

s. 

386.03 

MAIN  GEAR  C.G. 

TO 

FWO  C.G. 

M 

= 

41,70 

MAIN  gear  C.G. 

TO 

AVE  C.G. 

AVM 

29,20 

GROUND  LEVEL 

TO 

FWO  C.G. 

U 

= 

150.00 

GROUND  LEVEL 

TO 

ENGINE 

UE 

153,00 

BETWEEN  C.G. »S 

OF 

LEFT  K RIGHT  MAIN 

E 

= 

172,00 

TWIN  SPACING 

OF 

NOSE  TIRES 

SN 

= 

22.00 

TWIN  SPACING 

OF 

MAIN  TIRES 

SM 

0 .00 

TANDEM  SPACING 

OF 

MAIN  TIRES 

SNM 

= 

60.00 

SOIL  parameters 

CO^E  INDEX  CI=  192.00 

SOIL  TYPE  NTYP=  1 


ENGINE  PARAMETER 

THRUST  niEFFRENTIAL  PE= 


1.00 


MAX  STEERING 


tiO.OOMAX  SLIP  = 


0.00 


Figure  23.  Input  Data  for  Sample  Calculation. 
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In  the  original  program,  specific  values  of  steering  and  slip  were 
entered  as  half  second  functions  of  time,  with  many  runs  required,  and  since 
specific  values  were  not  really  known,  it  was  decided  to  use  40  time  pairs 
which  would  have  fractional  values  up  to  unity.  In  this  manner  a waveform 
could  be  entered  and  the  maximum  value  entered  with  one  card.  Thereafter 
all  inputs  could  have  the  same  waveform,  and  changes  could  oe  made  with 
only  one  card. 

Sample  output  data  for  the  inputs  of  Figure  23  are  shown  in  Figure  24. 
As  listed,  the  equilibrium  runwav  width  required  will  be  75.8  feet.  The 
results  shown  in  Figure  13  indicated  that  for  the  same  vehicle  with  60°  of 
nose  wheel  steering  the  field  width  would  be  106  feet.  This  difference  is 
not  due  entirely  to  the  subroutine  changes,  but  more  because  of  the  time  to 
enter  the  turn.  The  previous  run  was  for  a 10  second  entry  into  the  turn 
whereas  the  revised  data  are  for  a 2.5  second  entry. 

The  output  shown  reflects  several  changes  made  to  improve  the  output 
for  the  analyst.  The  pneumatic  trail  information  was  incorporated  by 
presenting  the  torque  generated  at  each  landing  gear.  For  example,  the 
turning  torque  listed  for  the  left  main  gear  reflects  the  total  torque  generated 
for  all  wheels  at  that  mount.  Because  of  the  use  of  lateral  force  ratio  in 
the  military  standard,  it  has  been  generated  for  each  gear  mount.  Lastly, 
because  of  the  need  to  recognize  the  soil  surface  damage  that  might  occur, 
instantaneous  sinkage,  and  rut  depth  are  computed. 

F.  TURNING  OPERATIONS  PARAMETRIC  STUDY 

The  digital  program  evolved  reflected  nearly  all  aspects  of  the  tire/ 
soil  interaction  that  could  be  analytically  incorporated.  At  this  stage  it  was 
desirable  to  use  the  program  and  compute  some  qualitative  results  as  a 
means  of  inferring  the  effects  of  the  many  parameters  that  influence  turning 
response.  The  number  of  parameters  is  considerable  when  we  consider, 
aircraft  type,  landing  gear  configuration,  tire  types,  soil  strength,  soil 
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TTME=  7,00  Xir  293»2  Y=  52.7  THRUST=  20345,3 

NOSE  GEAR  DATA 

TURN  ANGLE=  .21519  STEERING  ANGLE=  1.04716 
ORAG-rORCE"  = 1538.0  SICE  FORCE  = 6675.6 

DRAG  FORCE  RATIO=  .084  SIDE  RATIO  = .366 

TURNING  TORQUE  =12548.6  LATERAL  RATIO=  .376 

?=  .621  Z TURNEO=  1.034  RUT=  .574 

LEFT  MAIN  GEAR  DATA 

TURN  ANGLES  .04905  BRAKED  SLIP  PERCENTS  0.00 
- DRAG-  FORCE  = 4742.0  SIDE  FORCE  = 4557.5 

A 

DRAG  FORCE  RATIO=  .098  SIOE  RATIO  = .094 

TURNING  TORQUE  =31260.8  LATERAL  RATIO=  .136 

7=  1.125^-7  BRAKEO=  2.359  RUT=  1.310 

RIGHT  MAIN  GEAR  DATA 
TURN  ANGLES  .02935 

DRAG^  FORCE  = 9554.0  SIDE  FORCE  = 3755.0 

DRAG  FORCE  RATIO=  ,144  SIOE  RATIO  = .057 

TURNING  TORQUE  =11841.8  LATERAL  RATTO=  .155 

7=  ' 1.854  "7  TURNEO=“  2,050  RUT=  1,139 

AIRCRAFT  SIOE  FORCE  RATIO  = .077 

SUM  OF  FORCESs  -.00  SUM  OF  MOMENTS=  .00 

CENTER  OF  CURVATURE  AT  X=  12.57  Y=  340.97  RADIUS  342.40 

taur«  IBRIUM  REACHED, RUNNA-f  WIDTH  75,8 

Figure  24.  Output  Data  for  Sample  Calculations. 
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type,  gross  weight,  speed,  thrust  and  braking;  to  mention  a few  of  the  input 
parameters.  Similarly,  the  number  of  possible  outputs  are  considerable. 
Consequently,  it  was  decided  to  examine  specific  types  of  vehicles  for  fixed 
configurations.  The  three  types  selected  were  an  attack  aircraft,  a C-130 
configuration,  and  an  advanced  transport  configuration.  The  dimensional 
data  used  are  shown  in  Appendix  A.,  These  were  selected  because  the  attack 
vehicle  has  a smaller  gross  weight  with  all  single  wheel  landing  gear,  the 
C-130  has  a middle  gross  weight  with  twin  nose  gear  and  tandem  main  gear; 
and  the  advanced  transport  has  a greater  weight  with  twin  and  tcindem  main 
mounts. 

The  time  restraints  of  the  research  program  requireu  careful  selection 
of  the  number  of  parameters  and  parameter  variations  to  be  studied.  It  was 
decided  that  each  vehicle  would  be  examined  for  a fixed  gross  weight,  one 
taxi  speed,  one  tire  deflection,  one  soil  type  and  strength.  This  permits 
then,  evaluation  of  the  effects  of  steering  angle,  thrust  differential  and 
braking.  Steering  angles  of  20,  40,  cind  60  degrees  were  selected  to  provide 
data  indicative  of  current  nose  wheel  steering  limits.  Thrust  differentials 
of  unity  and  one -half  were  used.  Inboard  braking  of  0,  15,  30,  and  45  per- 
cent slip  were  used  to  generate  drag  ratios  up  to  0.  5. 

Several  outputs  are  possible  from  the  tabulated  data.  Those  most 
applicable  to  current  specifications  would  be  side  force  ratio  on  any  tire 
and  lateral  force  ratio  on  any  braked  tire.  Other  values  thought  to  be  critical 
were  runway  width  and  thrust  required  since  these  are  measured  of  how  well 
the  vehicle  performs  cr  can  perform.  These  data  are  all  plotted  on  one 
figure  to  show  the  effects  of  selected  inputs  on  desired  outputs. 

The  results  for  the  three  vehicles  are  shown  in  Figures  25  through 
33. 

All  vehicles  have  similar  performance  in  terms  of  the  apparent  curva- 
ture of  the  presented  lines  of  constant  parameter  values.  However,  each 
covers  a different  range  of  runway  width  and  force  ratio  magnitude.  The 


37 


Force  Ratios  {S/ 


/■ 


Maximum  Thrust  Available 


Figure  25.  C-130  Configuration  Performance  in  Cohesive  Soil  of 
CBR  6 at  5 Knots  and  Gross  Weight  of  133,000  Pounds 
with  Symmetrical  Thrust. 


Force  Ratios  (S/P  and  1>/P)  Thrust  (1000  Pounds) 


Figure  26.  C-130  Configuration  Performance  in  Cohesive  Soil  of 
CBR  6 at  5 Knots  and  Gross  Weight  of  133,000  Pounds 
with  50  Percent  Thrust  Differential. 
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Force  Ratios  (S/! 


Figure  28.  Advanced  Transport  Configuration  Performance  in  Cohesive  Soil  of  CBR  6 
at  5 Knots  and  Gross  Weight  of  216,000  Pounds  with  Symmetrical  Thrust. 


Figure  30.  Advanced  Transport  Configuration  Performance  Composite  for  Cohesive 
Soil  of  CBR  6 at  5 Knots  and  Gross  Weight  of  2l6,  000  Pounds. 


Force  Ratios  (S/P  and  L/P) 


Figure  31.  Attack  Vehicle  Performance  in  Cohesive  Soil  of  CBR  6 at  5 
Knots  and  Gross  Weight  of  31, 000  Pounds  with  Symmetrical 
Thrust. 
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Figure  33.  Attack  Vehicle  Performance  Composite  for  Cohesive 
Soil  of  CBR  6 at  5 Knots  and  Gross  Weight  of  31,  000 
Pounds. 


r.ttack  vehicle  data  presented  are  for  20,  30,  and  40  degree  nose  wheel 
steering  angles  because  the  60  degree  variation  could  not  be  calculated. 

The  program  prints  out  diagnostics  that  indicate  when  the  load  strength 
ratios  exceed  predictive  equation  limits,  turn  angles  are  too  great,  or  the 
iteration  cannot  converge  to  a solution.  For  the  attack  vehicle  all  occurred 
at  60  degree  nose  wheel  steering.  These  types  of  diagnostics  were  not 
included  in  the  original  program  and  it  was  then  necessary  to  look  at  the 
specific  values  printed  out.  A similar  response  was  calculated  for  the 
advanced  transport  configuration.  Consequently,  it  was  also  examined  for 
20,  30,  and  40  degree  nose  wheel  steering.  One  known  reason  for  the 
diffiulty  in  computing  large  steering  angle  response  is  that  the  turn  angle 
enters  into  the  calculation  of  turned  sinkage  as  indicated  by  measured  data. 
This  is  then  used  in  the  calculation  of  braked  drag  for  the  inboard  gear  and 
is  reflected  in  the  large  drag  force  generated  at  relatively  low  slip  values. 

The  relation  between  turned  sinkage  and  braked  slip  is  neither  known  nor 
measured.  However,  from  previous  works  it  was  assumed  that  braked 
response  would  be  calculated  from  instantaneous  sinkage  regardless  of 
whether  due  to  straight  ahead  or  turned  rolling. 

The  C-130  configuration  did  permit  calculations  over  the  desired 
range  of  pareimeter  values.  Whether  this  csui  be  attributed  to  wheel  base, 
track,  weight  or  tire  configuration  is  not  really  known  at  present.  There 
are  just  too  many  variables  that  influence  the  repponse  that  have  not  been 
examined.  The  C-130  plots  indicate  several  aspects  which  are  in  general 
true  for  the  other  two.  Starting  from  a nose  wheel  steering  only  configuration, 
inboard  braking  naturally  decreases  the  runway  width  required.  However, 
the  change  of  45  percent  slip  is  not  as  significant  as  that  gained  by  a 10 
degree  increase  in  steering  angle.  However,  the  10  degree  change  causes  a 
more  significant  change  than  a 50  percent  thrust  differential.  Since  the  main 
gear  lateral  force  ratio  is  primarily  dictated  by  the  braking,  it  appears  that 
the  best  method  of  achieving  reduced  structural  loads  into  the  gear  is  that 
of  not  using  the  brakes  unless  necessary. 
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The  one  aspect  that  is  difficult  to  generalize  is  that  of  the  nose  wheel 
side  force  ratio.  The  curves  that  are  plotted  above  the  abscissa  are  positive 
side  forces  indicating  that  the  nose  wheel  helps  to  turn  the  vehicle.  Those 
below  the  abscissa  indicate  that  the  nose  wheel  steering  angle  is  insufficient 
for  the  equilibrium  turn  generated.  The  nose  wheel  that  generates  a negative 
side  force  ratio  is  being  plowed  through  the  soil.  Therefore,  the  plot  for  a 
50  percent  thrust  differential  indicates  that  there  is  only  one  point  where  the 
steering  is  sufficient,  at  60  degrees  with  no  braking.  All  other  points  have 
insufficient  turning  for  a 50  percent  thrust  differential. 

The  thirust  curves  indicate  that  differential  thrust  not  only  decreases 
the  runway  widtli  required  but  does  so  at  reduced  thrust  required.  The 
danger  in  going  too  far  with  that  relation  is  that  as  the  thrust  required  goes 
down,  so  toe  does  the  available  thrust.  The  C-130  has  roughly  40,000  pounds 
of  thrust  available.  At  a minimum  radius  turn  the  thrust  required  is  Z4,  000 
pounds.  However,  in  order  to  get  a differential  of  50  percent,  the  outboard 
engines  must  develop  two-thirds  of  that  required,  16,  000  pounds  which  is  80 
percent  of  that  available  from  the  two  engines.  If  no  thrust  differential  exists 
for  the  same  steering  and  braking,  the  thrust  required  is  about  30,000  pounds 
or  75  percent  of  the  thrust  available  on  both  sides. 

The  plots  for  both  the  C-130  and  the  advanced  transport  indicate  that 
turning  is  best  accomplished,  in  terms  of  lateral  forces  on  the  main  gear, 
by  using  nose  wheel  steering,  then  differential  thrust,  and  lastly  braking. 

This  is  particularly  true  for  the  advanced  transport  configuration  shown  where 
a thrust  differential  essentially  treuislates  the  lateral  forces  ratio  curves 
without  increasing  them.  Therefore,  the  paragraph  of  MIL-A-00862A  per- 
taining to  turning  has  the  means  of  turning  the  aircraft  listed  in  the  optimum 
manner  by  stating  "unsymmetrical  thrust  or  nose  gear  steering"  first,  and 
"differential  braking"  last. 


It  is  difficult  to  discuss  at  length  any  additional  parameter  variations 
since  comparison  between  thrust,  braking,  and  steering  response  can  be  seen 
from  the  plots  for  the  specific  values  of  input  used.  What  happens  when  the 
gross  weight  changes,  soil  type  and  cone  index  changes,  or  velocity  increases 
i'j  not  known  exactly  although  other  sample  turns  have  provided  some  insight 
into  what  would  happen  qualitatively.  This  suggests  also  the  fact  that  the 
results  shown  are  indeed  qualitative,  \lthough  the  drags,  sinkages,  side 
forces,  etc. , are  all  calculated  from  predictive  equations,  these  were 
developed  from  laboratory  tests  of  wheels  not  complete  vehicles.  Therefore, 
the  shift  due  to  a given  braking  or  tlurust  is  probably  reasonably  correct  for 
a percentage  change,  but  may  not  create  exactly  27.4  feet  of  runway  width 
change.  Similarly,  the  C-130  will  for  given  conditions  have  a smaller  field 
width  required,  but  not  necessarily  50.0  feet. 


G.  AIRCRAFT  TURNING  SUMMARY 


Turned  tire  test  results  have  been  incorporated  into  a digital  program 
containing  drag  and  sinkage  prediction  equations  previously  evolved.  The 
equations  reflected  free  roiling  and  braked  drag  as  well  as  multiple  wheel 
effects  for  aircraft  tires  operating  in  cohesive  and  cohesionless  soils.  The 
turned  tire  data  were  used  to  add  side  force  curves,  pneumatic  trail,  in-line 
drag  modifications,  slip  due  to  turn  angle,  siiikage  due  to  turn  angle,  and  rut 
depth  prediction.  The  equations  of  motion  for  a rigid  vehicle  turning  at  a 
constant  velocity  are  solved  for  given  conditions  of  braking,  nose  wheel 
steering  and  thrust  differential.  The  program  generates  performance  data 
in  terms  of  path  over  the  ground,  side  and  drag  force  ratios,  lateral  force 
ratios,  turn  angles,  torque  due  to  pneumatic  trail,  instantaneous  sinkage, 
braked  sinkage  and  rut  depths. 


Several  plots  for  representative  aircraft  have  been  presented  to 
indicate  the  nature  of  tlie  turning  phenomena  in  terms  of  runway  width  . J 

required,  force  ratios,  and  thrust  required.  Some  information  computed  I ^ 

i 1 
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was  consistent  with  the  current  turning  criteria,  Tr.at  is,  the  C-130 
representation  generated  performance  curves  for  the  range  of  steering  angles 
permitted  u/  the  pilot's  handbook.  However,  the  others  could  not  be  taken 
to  their  limit  steering  angle  because  the  computed  turn  angles,  load  strength 
ratios,  or  other  internal  parameters,  v/ore  outside  of  the  range  of  the  pre- 
dictive equations.  When  attem.pting  to  compute  minimum  runway  width, 
the  existing  test  data  is  insufficient  and  the  program  has  difficulty  in  iteratir.g 
to  convergent  solution.  Practically  this  infers  that  it  is  difficult  to  prettict 
what  a tire  does  at  a turn  angle  of  70  degrees,  or,  what  really  happens  when  a 
landing  gear  nearly  pivots  about  a point? 


SECTION  in 


MULTIPASS  TIRE /SOIL  INTERACTION  ANALYSIS 

A.  GENERAL  AND  DEFINITIONS 

Multipass  criteria  for  the  operation  of  aircraft  over  remote  soil 
runways  have  been  previsonly  established.  The  Military  Vehicle /Aircraft 
Rut  Depth  Method  has  been  proposed  and  published  criteria  and  procedures 
for  calculating  flotation  are  available  in  "Aircraft  Ground  Flotation 
Investigation,  Part  I,  "AFFDL-TR-66-43  or  AFSCM  80-1  "Handbook  of 
Instructions  for  Aircraft  Design.  " 

The  current  research  program  was  designed  to  determine  whether  or 
not  new  criteria  could  be  developed  from  multipass  teat  data  which  would  be 
an  extension  of  the  previously  developed  predictive  equations.  Equations  and 
nomograms  are  available  to  calculate  instantaneous  sinkage  and  drag  for 
a rolling  single  wheel.  This  can  be  adjusted  to  reflect  twin  and  tandem 
configurations  as  well  as  braked  effects.  The  next  logical  steps  were  then 
to  determine  how  multiple  passes  of  the  rolling  or  braked  tire  over  alternating 
but  adjacent  paths  could  influence  sinkage,  drag  and  rut  depth.  If  trends 
could  be  established  as  functions  of  the  number  of  the  pass,  v/here  it  occurred, 
and  whether  or  not  it  were  braked;  relative  to  first  pass  calculated  response; 
then  predictive  equations  would  be  available  to  better  establish  "failure 
criteria". 

Current  multipass  criteria  are  functions  of  drag  developed  and 
cumulative  rut  depth.  WiP  the  drag  exceed  the  capability  of  the  vehicle  to 
maneuver  over  the  soil?  Will  the  cumulative  rut  depth  cause  excessive 
drag  or  interference  between  vehicle  structure  and  the  soil  surface?  If 
the  data  permits  reasonable  prediction  of  cumulative  rut  depth  and  drag  as 
a function  of  the  multipass  history  (rolling,  braked,  same  lane, adjacent 
lane),  then  it  would  be  possible  to  better  relate  the  physical  characteristics 
of  the  vehicle,  its  mission,  cind  flotation  criteria. 
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In  the  following  sections  it  is  well  to  remember  that  the  predictive 
equations  that  are  available  are  based  upon  instantaneous  sinkage  calculated 
for  the  tire  type,  applied  load,  cone  Index  and  soil  type.  Measured  test 
data  does  not  yield  instantaneous  sinkage  and  hence  it  is  necessary  to  relate 
it  to  another  parameter  such  as  axle  movement  for  the  AEWES  data.  After 
instantaneous  sinkage  is  found  to  relate  to  measured  drag,  then  it  is  ailso 
correlated  with  rut  depth,  another  measured  parameter.  Therefore,  in  the 
analysis  of  the  data  it  will  be  shown  that  indeed  sinkage,  axle  movement, 
and  rut  depth  are  different  quantities  which  ctm  be  related  not  only  for  first 
pass,  but  multipass  operations. 


B.  TEST  PROGRAM  AT  AEWES 


Test  Matrix 

A series  of  tests  were  conducted  at  the  Army  Engineers  Waterways 
Experiment  Station  (AEWES)  under  the  supervision  of  University  of  Dayton 
personnel.  It  was  desirable  to  have  three  basic  types  of  tests.  These  were 
free  rolling  sinlge  wheel  over  one  lane,  rolling  and  bradced  over  alternating 
adjacent  lanes,  and  twin  tires  free  rolling  over  alternating  adjacent  lanes. 
The  tests  were  to  be  conducted  for  at  least  two  types  of  tires,  two  cone 
index  values  and  two  applied  loads.  The  original  test  outline  is  shown  in 
Figure  34.  The  pass  distribution  for  alternating  paths  are  shown  in 
Figure  35. 


From  these  figures  a test  matrix  was  evolved  to  specify  the  testing 
to  be  accomplished.  The  first  page  of  the  matrix  is  shown  in  Figure  36. 
The  variable  of  greatest  importance  was  soil  type.  That  is,  the  variable 
which  could  be  changed  least  frequently.  Hence,  the  test  matrix  was  set 
up  to  run  all  high  cone  index  soil  * sts  first.  Secondly,  single  tire  tests  of 
one  tire  type  precedes  a second  tire  type,  with  variable  applied  loads. 
Therefore,  the  test  plan  was  to  progress  with  minimal  cheinges  from  one 
test  to  the  next  and  one  series  to  the  next. 
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Figure  34.  Multipass  Tire  Test  Outline. 


Passes  Distribution  for  Test  Part  2,  Twin  Tire  Configuration, 
for  18  Passes  as  Shown. 
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Figure  35.  Passes  Distribution  for  Test  Part  2 and  Test  Part  3. 
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MULTIPASS  TIRE  TEST  PROGRAM  TEST  MATRIX 


Test 

No. 

Tire 

Designation 

Cl 

Avg 

Vertical 

Load 

lbs 

Number 
of  'fires 

Pass  Configuration 
Path  Draking 

Rut  Depth 
and 

Cone  Index 

1 

7:00-6.6  PR 

IB 

800  to  900 

Single 

All 

Free 

Yes 

■ 

B 

1 

in 

Same 

Rut 

Rolling 

Yes 

B 

B 

800  t 

0 900 

Yes 

11 

■mm 

l^mm 

■ 

1100  to  1200 

Yes 

m 

B 

B 

1 

1 

1 

'{c.s 

16 

17 

18 

19 

20 

1 

1 

1 

1 

1 

1 

All 

in 

Same 

Rut 

Free 

Rolling 

Yes 

21 

■ 

1 

R 

Yes 

22 

I 

*» 

B 

23 

1 

3 

B 

24 

■' 

1 

B 

25 

m 

1 

B 

Yes 

26 

1 

1 

R 

27 

■ 

2 

R 

28 

1 

2 

B 

29 

m 

1 

R 

30 

■ 

1 

R 

Yes 

B 

3 

B 

■ 

B 

3 

R 

B 

3 

R 

■ 

B 

2 

R 

I9BI 

m 

B 

1 

1 

R 

Yes  « 

2^ 

1100  to  1200  1 

SintJle  1 

1 

3 1 

Notes: 

Cl  of  55  - CBR  of  1.7 


Sec  Figures 
1 and  2 for 
Code  to 
Particular 
Paths 


♦ Measured 
after  the 
Numbered 
Test  or  When 


Failure 


The  number  of  passes  over  a particular  path  had  been  selected  as 
indicative  of  the  distribution  that  would  exist  about  the  centerline.  The 
number  of  rolling  versus  braiked  passes  was  selected  to  provide  sufficient 
data  indicative  of  the  acceleration  of  "failure"  due  to  both  rolling  and 
braking  action.  Once  these  were  established,  the  selection  of  alternating 
rolling  or  braked,  this  path  or  that,  were  determined  in  a purely  random 
selection  process. 

Failure  Criteria 

Prior  to  any  testing  it  was  necessary  to  establish  failure  criteria 
guidelines  for  the  tests.  Two  criteria  were  established.  Instantaneous  sinkage 
divided  by  the  tire  diameter  was  restricted  to  0.06.  This  was  established  by 
examining  the  drag  ratio- sinkage  ratio  of  a single  wheel  on  soil  shown  in 
Figure  37.  Assuming  that  most  cargo  aircraft  have  a thrust  to  weight  ratio 
of  0.25  to  0.  30,  takeoff  would  be  impossible  if  the  R/P  ratio  were  greater 
than  approximately  0.25.  This  establishes  0.06  as  the  limiting  value  for  Z/D. 
Hence,  for  a 8:50  tire  the  instantaneous  sinkage  would  be  approximately  one 
and  one -half  inch. 

Cumulative  rut  depth  to  tire  diameter  ratio,  Z^^/D  was  restricted  to 
0.20.  This  approximates  the  value  of  one  half  to  two  thirds  of  the  free 
height  of  the  aircraft  tire,  or  approximately  five  inches  for  the  larger  tire. 

Test  Measurements 

Tests  were  conducted  using  Type  III  aircraft  tires;  7:00-6,  6PR  and 
8:50-10,  8PR,  both  operating  at  deflections  of  35%.  All  tests  were  conducted 
at  forward  speeds  of  10  feet  per  second.  The  following  measurements  were 
taken: 
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Figure  37.  Drag  Ratio -Sinkage  Ratio,  Single  Wheels  on  Soil. 


Vertical  Load  (each  pass) 

Drag  Load  (each  pass) 

Axle  Movement  (each  pass) 

Rut  Depth  (each  pass) 

Rut  Profile  (after  test,  5th,  10th,  and  either  last 
or  failure  pass) 

Velocity  (each  pass) 

All  data  collected  were  reported  in  Reference  18.  The  following  analyses 
select  typical  test  results  from  that  report. 


C.  SINGLE  WHEEL  FREE  ROLLING  TESTS 


Test  Part  1 of  the  test  matrix  called  for  four  series  of  tests  to  be 
conducted  using  two  tire  types,  two  soil  strengths,  two  load  levels  and  ten 
passes  in  the  same  rut.  Each  test  was  exa.mined  to  determine  how  well  the 
first  pass  data  matched  existing  techniques,  and  then  to  establish  a relation 
between  drag  and  rut  depth  as  a function  of  the  pass  number  and  first  pass 
calculated  valves. 

Test  series  A74-0001-3  was  conducted  using  a 7:00-6  tire  at  35% 
deflection.  The  raw  data  are  plotted  as  shown  in  Figure  38,  and  as  tabulated 
below. 


TEAt. 

74. 

TEST  1. 

SOIL 

3. 

fASS 

SO. 

OESI&ED 

SLIP 

Z 

ACTUAL 

SLIP 

Z 

LOAD 

LBS 

DRAG 

LBS 

AXLE 

MOVIMENT 

m. 

AXLE 
VELCCITY 
FT/ SEC 

TORQUE 

FT-LB 

1. 

1.2 

798.3 

-64.7 

0.48 

10.219 

0.0 

2. 

2.2 

747.1 

-66,3 

0.34 

10.208 

0.0 

3. 

1.7 

754.0 

-66.0 

0.55 

10.240 

0.0 

A. 

1.4 

753.2 

-70.0 

0,59 

10.230 

0.0 

5. 

1.8 

738.7 

-70.6 

0.71 

10.298 

0.0 

6. 

2.0 

762.0 

-7:..3 

0.79 

10.185 

0.0 

7. 

2.0 

753.4 

-7fc.5 

0.86 

10.301 

0.0 

8. 

1.2 

759.8 

-77.9 

0.92 

10.346 

O.C 

9. 

2.2 

753.5 

-78.2 

0.99 

10.354 

0.0 

.0. 

2.2 

764.3 

-77.2 

1.08 

10.408 

0.0 

58 


The  cone  index  values  were  measured  before  the  test  began  and  after 
the  first,  fifth,  and  tenth  pass. 

The  first  pass  generated  an  eixle  movement  that  appears  to  be  above 
the  general  trend  of  the  last  eight  passes.  We  would  anticipate  that  first 
pass  effects  would  differ  from  later  pass  effects,  but  not  as  drastically 
different  as  shown.  Consequently,  the  first  pass  data  was  closely  examined 
to  determine  how  well  it  could  be  explained  and  understood. 

For  the  7:00  tires,  the  rigid  surface  contact  area  is  calculated  to  be 
43.4  square  inches.  Using  the  nomogram  shown  in  Figure  39,  the  sinkage 
and  drag  were  calculated  for  the  first  pass.  An  applied  load  of  798  pounds 
with  a cone  index  of  43.  0 generates  a sinkage  of  0.  75  inches  and  an  R/P  of 
0.  157.  Hence  calculated  drag  would  be  125  pounds  as  compared  with  65 
measured. 

The  rut  profile  data  indicates  that  the  rut  depth  measured  after  the 
first  pass  was  0.097  (average  for  three  locations).  The  axle  movement 
measurements  show  a motion  of  0.48  inches.  From  AEWES  data  relating 
sinkage  to  axle  movement  (Figure  40),  the  relation 

Z =2/3  Z 
a 

was  used  for  35%  tire  deflection.  The  calculated  axle  motion  would  then 
have  been  0.  50  inches. 

The  relation  between  instantaneous  sinkage,  Z , and  rut  depth  comes 
from  Figure  41,  data  collected  at  the  University  of  Dayton.  The  curve  shown 
is: 


Z = (2.4  Zj^) 


1/2 


The  calculated  rut  depth  would  then  be,  Z = 0.21  inches,  therefore; 

Xv 
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FLOTATION  INDEX  - OPERATIONS  INDEX  NOMOGRAPH 
SINGLE  WHEEL 


Figure  39.  Nomograph  for  Calculation  of  Instantaneous  Sinkage  and  Drag. 


AEWES  Data  Relating  Sinlcage  to  Axle  Movement, 


Figure  41.  Relation  Between  Instantaneous  Sinkage,  Z,  and  Rut  Depth  -UD  Data 
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Measured  Data 

Calculated  Data 

z 

(None) 

0.  75 

^A 

0,48 

0.  50 

^R 

0.097 

0.  21 

R 

P 

0.080 

0.  157 

The  sinkage  calculated  does  agree  with  th 
ates  drag  that  is  too  high  and  a rut  depth  that  is 

e axle  movement, 
too  high. 

but  gener- 

The  "linear"  portion  of  the  axle  movement  occurs  over  the  range  of 
3 to  10  passes.  In  this  area  the  average  drag  is  73.9  pounds  and  average 
applied  load  is  758,6  pounds.  The  nomographic  solutions  for  drag  and 
sinkage  using  a P/CI  of  16.  6 are; 

Z =0.^3 

R/P  = 0.095 

Using  a similar  approach  as  before,  the  calculated  axle  movement  would  be 
0.29  inches  and  rut  depth  would  be  0,0775, 

Therefore; 

Measured  Data 

Calculated  Data 

Z 

(None) 

0.43 

0.48 

0.29 

0.097 

0.0775 

R 

P 

0.080 

0.095 
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The  measured  data  are  those  collected  for  a firs*  pass  which  had  an 
applied  vertical  force  greater  than  desired.  The  niean  test  value  for  vertical 
force  was  about  760  pounds  but  the  first  pass  was  conducted  at  800  pounds. 
The  calculated  column  indicates  what  would  have  happened  if  the  "correct" 
mean  value  hau  been  used  for  the  first  pass. 

The  calculated  sinkage  cannot  be  validated  with  measured  data  because 
there  is  none.  The  axle  motion  similarly  cannot  be  compared  with  but  does 
indicate  that  theoretically  the  measurement  would  have  been  less  than  that 
observed  and  more  nearly  that  measurement  on  the  second  pass.  The  calcu- 
lated rut  depth  0.080  is  more  closely  related  to  the  measured  0,097,  and 
calculated  drag  would  be  72. 1 pounds  for  a very  close  correlation. 

Therefore,  the  first  sample  of  data  indicate  that  first  pass  calculations 
of  rut  depth  and  drag  do  approximate  those  measured  and  that  the  drag  over 
the  ten  passes  does  not  change  appreciably.  Similar  procedures  were  used 
for  all  four  single  wheel  multipass  tests  and  the  results  are  tabulated 
below. 


Test  Series 
Number 

Rut  Depths 

Mean 
T .nads 

P 

1 ■ 1 ■ i.i 

First  Pass  Calculated 

^A 

Z 

D 

1st 

5 th 

10th 

P 

"d 

Cl 

A74-0001-3 

0.  10 

0.  52 

0.98 

760 

74 

16.6 

0.  52 

0.78 

0.08 

84 

A74-0002-3 

0.  31 

1.41 

2.  36 

1060 

160 

22.8 

0.75 

0.50 

0.23 

159 

A74-0003-3* 

0.  !1 

0,  34 

0.60 

1170 

82 

22.4 

0.41 

0.27 

0.07 

82 

A74-0010-3* 

0.  16 

0.  57 

1.22 

1050 

105 

31.4 

0.73 

0.48 

0.22 

116 

* 8:50-10 

"ires 
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In  all  tests  it  appears  that  if  first  pass  calculation  procedures  are  used 
to  calculate  the  drag  and  rut  depth,  that  the  computed  and  mtpsured  will  be 
within  reasonable  agreement. 

The  second  portion  of  tlie  analysis  was  to  determine  whether  or  not 
predictive  techniques  could  be  evolved  to  calculate  rut  depth  as  a function 
of  number  of  passes. 

Cumulative  rut  depth  as  a function  of  number  of  passes  is  plotted  on 
the  figure  on  the  following  page-.  The  ordinate  reflects  the  ratio  of  measured 
rut  depth  to  first  pass  calculated  rut  depth.  This  was  done  since  we  are 
after  relations  which  can  be  used  to  predict  based  upon  first  pass  calculated 
values.  The  curve  indicates  that  there  is  a spread  in  the  data,  and  that  the 
rut  ratio  is  highest  for  the  smaller  tire.  For  each  tire  the  ratios  decrease 
with  increasing  P/CI  values. 

Various  parameter  relations  were  examined  to  determine  a commonality 
between  all  tests  which  could  be  used  to  evolve  a predictive  equation.  The 
result  can  be  expressed  as: 


'RIC 


^ =16.5 


where:  is  the  10th  pass  rut  depth 

Z is  the  first  pas^  calculated  rut  depth 

R.XC> 

P is  the  vertical  load  to  cone  index  ratio 
Cl 

A is  the  ratio  of  footprint  area  to  cone  in 
Cl 


Z - Z 


r J is  the  rebound  over  instantaneous  sinkage  for  the 
calculated  first  pass 
1C 

^ is  the  inverse  of  footprint  length. 

A/ 
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PASS  NUMBER 


Figure  42.  Variation  of  Cumulative  Rut  Depth  Z divided  by  First 
Pass  Calculated  Rut  Depth  Z 

^CI 


2 


O 


7:00  Tire  P/CI  = l6.6 
7:00  Tire  P/CI  = 22.8 
8:50  Tire  P/CI  = 22.4 
8:50  Tire  P/Cl  = 31.4 


Using  these  relations  for  the  mean  tabulated  values  of  P/CI  and 
appropriate  values  of  A and  I , the  computed  10  pass  rut  depths  are: 


Test  Series 
Number 

Rut  Depths 

Measured 

Calculated 

A74-0001-3 

0.98 

0.97 

A74-0002-3 

2.  36 

2.46 

A74-0003-3 

0.60 

0.63 

A74-0010-3 

1.  22 

1.  13 

If  it  is  at  sumed  that  rut  depth  increases  linearly  with  number  of 
passes,  then  th«’  rut  depth  can  be  calculated  from  the  expression: 

2n  (n-i) 

where  N = the  number  of  the  pass 

is  the  first  pass  rut  depth 

ZjP  is  computed  from  the  above  equation. 

A more  rigorous  analysis  could  be  conducted  to  better  fit  the  data,  but 
it  is  believed  that  it  is  not  warranted  at  the  present  time.  The  scatter  in 
the  rut  depth  data  as  well  as  the  differences  between  measured  axle  move- 
ment  and  computed,  drag  observed  and  computed,  variation  in  cone  index 
readings,  etc. , all  indicate  that  for  only  ten  passes,  the  relation  found 
should  be  sufficiently  accurate.  A plot  of  the  calculated  data  are  shown  on 
the  following  page  in  comparison  with  measured. 

The  results  of  the  analysis  indicate  that  first  pass  sinkage  calculations 
cein  be  used  to  compute  the  rolling  drag  and  cumulative  rut  depth  for  rolling 
tires.  The  sinkage  and  drag  calculations  are  from  previous  relations  valid 
for  many  tire  types  and  sizes.  The  cumulative  rut  data  is  based  upon  7:00-6 
and  8:50-10  tires  at  35%  deflection  operation  over  clay  with  cone  index 
values  of  37  to  55. 
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SINGLE  WHEEL,  RANDOM  BRAKING  AND  PATH 


Test  part  3 of  the  test  matrix  called  for  four  tests  to  be  conducted 
with  two  tire  types,  soil  strengths  and  load  levels.  Sixteen  passes  were 
made  randomly  alternating  over  three  adjacent  paths  with  or  without 
braking. 


Preliminary  Analysis 


Tests  identified  as  series  A74-0009-3  were  first  examined.  The 
preliminary  information  is  shown  below. 


YiriR 

74. 

TEST  9, 

SO 

DESIRED 

ACTUAL 

PASS 

SLIP 

SLIP 

LOAD 

HO., 

% 

LBS 

1. 

0.  0 

2.  3 

1 176.  4 

2. 

35.  G 

34.  0 

11  S3, 2 

3„ 

35.  0 

33.  9 

1 176.  4 

4., 

35.  0 

34..  7 

1 174.  9 

5.. 

0.  0 

2.  1 

1163.  2 

6.. 

0.  0 

2.  t. 

1171.3 

7. 

0 0 

5.  4 

11S4.  4 

8. 

35.  0 

34.  S 

11S5.  4 

9. 

0.  0 

2.  S 

1164.  5 

10., 

0.  0 

2.  6 

1164.7 

11. 

0.  0 

3.  7 

1172.  7 

12. 

35.  0 

~ ' 0 

1173.  7 

13. 

35.  C 

" r. 

1154  0 

14. 

0.:  i 

4.3 

1 1 64 , 0 

15. 

0.  0 

3.  7 

1 1 37 . 6 

16. 

35.  0 

35..  2 

1170.  7 

AXLE 

AXLE 

DRAG 

MOVE H EH T 

VELOCITY 

TOROUE 

LBS 

IH, 

FT/SEC 

FT-LB 

TOO.  3 

0.  21 

10. 357 

0.  0 

505.  5 

0.18 

10. 351 

-340  5 

516.  4 

0,44 

10,. 333 

-348.  3 

522.  7 

0.13 

10,449 

-338.  4 

112.4 

0.  60 

10,347 

0.  0 

106.  S 

0.  74 

10. 3S5 

0.  0 

157.  7 

-0,02 

10  350 

0.  0 

566.  2 

0.31 

10. 513 

-353.  2 

122.  0 

0.  57 

10,425 

0.  0 

122,  1 

0.  SO 

1 0. 377 

0.  0 

143,6 

0.  16 

10. 369 

0.  0 

607,3 

0.  77 

10.  466 

-388.  4 

5S0 . 4 

1.22 

10  52 S 

-378.  7 

164  7 

0 05 

10. 547 

0.  0 

142,9 

1 20 

10. 512 

0.  0 

5 £"5.  2 

0.  23 

10  56-! 

-328.  6 
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This  was  then  retabulated  in  the  form  shown  below. 


Axle 

Pass  Braked  or  Movement 


Number 

Path 

Rolling 

Path  1 

Path  2 

Path  3 

1 

1 

R 

0.90 

- 

- 

2 

2 

B 

- 

0.85 

- 

3 

3 

B 

- 

- 

0.81 

4 

1 

B 

0.84 

- 

- 

5 

1 

B 

1.43 

- 

- 

6 

1 

R 

1.98 

- 

- 

7 

2 

R 

- 

0.78 

- 

8 

2 

B 

- 

1.43 

- 

9 

1 

R 

2.06 

- 

- 

10 

1 

R 

2.  38 

- 

- 

11 

3 

B 

- 

- 

0.  72 

12 

3 

B 

- 

- 

1.45 

13 

3 

R 

- 

- 

2.  32 

14 

2 

R 

- 

1.76 

- 

15 

1 

R 

1.45 

- 

- 

16 

1 

B 

2.05 

- 

It  is  interesting  to  note  that  the  axle  movement  does  not  monotonically 
increase  along  any  one  path.  Apparently  the  axle  movement  over  one  path 
increases  as  the  wheel  passes  directly  over  the  same  rut,  but  decreases  if 
there  are  intervening  passes  on  adjacent  lanes.  This  suggests  that  direct 
application  one  one  lane  causes  a heaving  on  the  adjacent  lane  which  could 
or  could  not  vary  with  the  free  rolling  or  braked  nature  of  the  pass. 
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The  specific  tests  series  mentioned  was  examined  in  detail  as  dis- 
cussed in  Appendix  B.  For  the  first  analysis  it  was  assumed  that  axle 
movement  would  be  indicative  of  both  instantaneous  sinkage  and  rut  depth. 
This  was  a gross  simplification  but  was  used  to  determine  whether  the 
motion  could  be  predicted  or  not,  could  the  variation  be  rationally  explained 
or  would  this  difference  in  measurements  be  within  measurements  and 
integration  error  bands. 

The  approach  taken  was  to  assume  that  the  axle  movement  of  one  path 

is  the  summation  of  the  effects  of  all  prior  passes.  Also  it  was  assumed  that 

only  adjacent  lanes  would  influence  one  another.  Tl._  .efore,  if  the  center 

lane  had  a direct  application  of  a rolling  wheel,  the  movement  would  be  as 

shown  below  where  A implies  applied  rolling  and  N near  rolling  effect. 

R R 

Paths 

12  3 Measurement 

"Sinkages”  A N N 0.90" 

R R R 

Since  both  lanes  1 euid  3 are  adjacent,  some  near  rolling  effects  would 
be  generated.  These  are  defined  as  The  second  pass  was  a braked  pass 

over  the  second  path  and  the  cumulative  motxcn  would  be  , 

Paths 


1 

2 

3 

Measurement 

Pass  1 

0,90" 

Pass  2 

"■b 

0.85" 

In  this  case  A implies  applied  braking,  N , near  brsiking  effect,  and  F 

JD  O Sj 

for  braiking.  Notice  that  the  measured  axle  movement  is  only  for  the  lane 
of  directly  applied  load. 
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If  there  are  near  and  far  fields  effects  for  rolling  and  braking  that  are 
constants,  (fixed  applied  load,  cone  index  and  slip),  then  it  would  theoretically 
be  possible  to  sum  the  total  sinkage  down  to  any  particular  pass  and  calculate 
what  the  measurement  should  be.  There  are  16  passes  made  and  six  unknowns. 
It  should  be  possible  to  evaluate  the  constants. 

This  process  was  followed  in  the  detailed  analysis.  Additional  work 
conducted  had  shown  that  far  field  effects  were  very  small  indeed  and  con- 
sequently only  near  field  effects  were  evaluated.  Starting  with  the  first  pass 
and  sequentially  proceeding  down  the  tests  generating  simultaneous  equations, 
it  was  found  that  there  weren't  any  constant  values.  It  was  found,  however, 
that  there  were  trends  which  could  be  estab'.is’ied  to  explain  the  data. 

For  the  7:00-6  tire  rolling  over  clay  with  a cone  index  of  approximately 
43,  making  passes  separated  by  one  tire  width,  the  following  could  be  used 
to  predict  axle  motion 

n m k 1 


where 


II 

0.89 

= 0.89 

^R 

-0.  06 

R2 

0.  55 

^B2 

= 0.  58 

^B 

= 

-0.47 

R3 

0.  32 

^B3 

= 0. 66 

The  coefficients  infer  that  there  is  some  predictable  pattern  but  that  the 
true  nature  is  not  well  defined.  It  appears  that  the  first  pass  creates  more 
sinkage  than  successive  passes,  and  that  successive  sinkages  may  or  may 
not  decrease  depending  upon  rolling  or  braking.  It  also  infers  that  rolling 
or  braking  on  an  adjacent  path  provides  a decrease  in  rut  depth.  Near  braking 
of  -0.47  indicates  that  a braked  pass  raises  the  rut  by  approximately  one 
half  of  the  directly  applied  rut  depth. 


Detailed  Analysis 


There  were  four  test  series  conducted  using  single  wheel  alternating 
paths  and  braking.  Each  test  had  available  cumulative  rut  depth  profile  data 
over  the  lane  just  tested.  Similarly  each  test  had  available  the  integrated 
motion  observed  for  the  same  test.  The  profile  data  was  collected  over  14 
stations  and  averaged  to  generate  an  average  rut  depth  These  data  were 
plotted  for  each  test  as  shown  in  Figures  44,  45,  46,  and  47. 

In  the  figures  each  pass  is  identified  as  to  whether  it  was  braked  or 
rolling.  Also,  the  plot  makes  clear  those  passes  made  specifically  over  the 
middle  path  which  is  the  path  of  interest.  The  first  thing  noticeable  is  that 
rut  depth  and  axle  movement  are  indeed  approximately  equal.  Secondly,  in 
order  for  the  rut  depths  to  vary  as  shown  it  quickly  becomes  apparent  that 
the  indirect  application  must  indeed  create  a heaving  as  previously  suspected. 
If  the  first  pass  is  direct  rolling  and  the  second  and  third  indirect  braking, 
then  how  can  the  cumulative  rut  depth  at  the  end  of  the  fourth  pass  be  approxi- 
mately equal  to  the  rut  depth  at  the  end  of  this  first  pass?  If  the  direct 
braking  causes  a rut  depth  of  3/4  inch,  as  indicated  by  the  first  and  fifth 
pass,  tlien  the  rut  depth  at  the  end  of  the  third  pass  would  have  to  be  near 
the  original  ground  level.  Hence  two  indirect  braked  passes  did  raise  the 
soil  the  equivalent  of  one  direct  rolling  pass. 

At  this  point  there  are  many  approaches  possible  but  all  theoretical 
possibilities  are  not  practical.  For  example,  it  is  possible  to  evaluate  the 
change  of  rut  depth  as  a function  of  number  of  the  pass  or  number  of  the 
braked  pass  orly,  or  rolling  pass  only,  and  generate  a change  in  rut  as  a 
function  of  several  parameters.  However,  the  amount  of  data  available 
doesn't  permit  this  type  of  rigorous  approach.  The  first  rolling  pass 
generates  0.9  inch  rut,  the  second  creates  a change  of  0.  3 inches,  and  the 
third,  0.  3 inches.  Was  this  an  initial  step  with  subsequent  rut  increases 
at  constant  values?  Or  did  the  accuracy  of  the  measurements  and  inte- 
grations generate  observations  of  a geometric  series  such  as  0.9  inch,  0.  54 
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Figure  46.  Cumulative  Rut  Depth  and  Axle  Motion  Versus  Number  of  Passes 
for  Center  Lane. 


Cumulative  Rut  Depth  and  Axle  Motion  Versus  Number  of  Passes 
for  Center  Lane. 


inch,  and  0.324  inch?  With  the  limited  data  and  referring  back  to  the  single 
wheel  free  rolling  data,  it  seemed  reasonable  to  assume  that  the  first  pass 
does  create  a unique  rut,  and  that  all  subsequent  changes  are  constants 
independent  of  the  number  of  the  pass.  If  a direct  rolling  rut  is  developed 
on  the  tenth  pass,  the  change  in  rut  depth  will  be  the  same  as  for  the  second 
pass.  The  braked  rut  increment  will  be  a different  value,  but  have  that  same 
value  for  all  braked  passes. 

For  the  data  of  Figure  44  the  coefficients  for  direct  and  indirect  rolling 
and  braking  are 


where 


*^1 

+0.25 

+0.00 

+0.75 

- 0.  35 

is  rolling 

direct 

is  rolling 

indirect 

is  braked 

indirect 

is  braked 

indirect. 

Using  these  numbers  the  predictive  cumulative  rut  deplh  for  lane  one 
would  have  been  as  shown  in  Figure  48.  Similar  procedures  were  used  to 
develop  coefficients  for  each  test  series.  The  values  are  shown  in  the 
following  section  where  they  are  used  in  the  predictive  auialysis. 

E.  MULTIPASS  PREDICTIVE  PROCEDURES 

It  was  originally  desired  to  have  a computer  routine  which  would 
perxnit  prediction  of  multipass  effects  based  upon  test  data.  Due  to  delays 
in  obtaining  the  data,  it  was  not  possible  to  generate  the  digital  routine, 
but  it  was  possible  to  evolve  the  predictive  equations  for  the  data  that  were 
available. 


IS. 
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Path  §\  on  Passes  1,  4,  5,  6,  9,  10,  15,  16;  Hence,  Available 
Rut  Depth  Versus  Predicted  Using  Coefficients  of  page  85, 


Previous  sections  indicated  the  results  of  tests  conducted  to  determine 


the  performance  of  a single  wheel  passing  over  a soil  surface  for:  (1) 
multipass  passes  over  the  same  path  while  free  rolling,  and  (2)  multiple 
passes  over  alternating  paths  while  randomly  braked.  The  data  discussed 
in  the  following  section  does  not  reflect  the  additional  data  collected  at  the 
University  of  Dayton  which  may  extend  the  applicability  of  the  predictive 
equations. 

Through  analysis  of  the  data,  it  has  been  shown  possible  to  predict 
cumulative  rut  depth  and  drag  as  a function  of  the  time  history  of  multipass 
operation.  If  the  tire  passes  over  the  same  path,  or  an  adjacent  path,  and 
rolls  or  is  braked,  it  is  possible  to  accumulate  the  effects  of  the  passes. 
The  steps  are  as  follows, 

1)  Calculate  first  pass  sinkage  and  drag  as  found  using  existing 
techniques. 

2)  Calculate  the  braked  drag  using  existing  techniques. 

3)  Calculate  first  pass  rut  depth. 

4)  Calculate  incremental  rut  depth  changes  for  rolling  passes. 

5)  Calculate  incremental  rut  depth  changes  due  to  direct  braking  or 
indirect  rolling  and  braking. 

6)  Define  history  of  multipasses. 

7)  Ceilculate  cumulative  effects. 


F.  DISCUSSION 

Calculation  of  First  Pass  Sinkage  and  Drag 

For  a given  tire  with  known  applied  load,  the  sinkage  and  drag  are 
calculated  by  the  equations: 

-f-  = 0. 11  +0.  33  (a/CI)  for 

X/ 

Y = -0.  03  +0. 19  (a  /CJ.)  for^ 

and  = 0.18  +3.23  (Z/D) 

The  equations  for  Z/t  are  specifically  for  cohesive  soil  only.  This  is 
consistent  with  all  other  data  discussed  in  this  report. 

The  calculations  can  be  accomplished  either  using  the  above  equations 
or  the  nomogram  which  is  available. 

Calculation  of  Braked  Drag 

Braked  drag  is  calculated  by  first  finding  the  rolling  tire  sinkage  aoid 
drag  using  the  equations  above.  The  ratio  of  fully  braked  sinkage  to  rolling 
sinkage  is  estimated  to  be  3.0.  Therefore, 


Z 


in  the  ranges  of  0.  6 to  0.  7. 
in  the  ranges  of  0.2  to  0.6 


where  S is  the  slip  desired. 


'■k 


The  values  of  and  S are  then  substituted  into 

B 


0.018  +3.23 


+ 0.09 


First  Pass  Rut  Depth 

Previous  AEWES  and  UD  teat  results  had  shown  that  it  was  possible 
to  relate  axle  motion,  instantaneous  sinkage  and  rut  depth  by  the  following 
expressions. 

Z = 2/3  Z Z < 3.0  inches 

A 

2 

and  Z =2.4Z  Z <2.5  inches 

rut  rut 

The  first  expression  is  a linear  fit  to  the  curve  of  Z^  versus  Z for 
a tire  having  35%  deflection  as  reported  in  AEWES  TR  No.  3-516  by  Smith 
and  Frietag.  The  second  expression  is  quadratic  fit  to  rebound  test  results 
measured  at  UD  with  limits  of  2.  5 inches  rut  depth. 

From  the  most  recent  test  program  it  was  found  that  the  curves  had 
to  be  modified  in  order  to  match  measured  aixle  motions  and  rut  depths. 
The  modifications  require  that; 


Z :=  0.55  + 1.09  Z^  Z > 1.1" 

Jbv 

Z^  = 2.4  Z„  Z < 1.1" 

K. 


Therefore,  If  first  pass  sinkage  is  known,  the  rut  depth  can  be  calculated. 
The  equation  was  found  using  the  axle  motion  equation 


Z*  =2/3  Z 
A 
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Lasts' 


\ 
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since  instantaneous  sinkage  is  not  a measured  parameter.  It  was  found  that 
computed  instantaneous  sinkage,  Z,  generates  axle  motions  and  rut  depths  which 
agree  with  those  measured  for  the  first  pass. 

Incremental  Rut  Depth  - Rolling 

The  free  rolling  multipass  data  indicated  that  it  was  possible  to  cal- 
culate cumulative  rut  depth  based  upon  first  pass  data.  Specifically, 


The  tenth  pass  rut  depth  is  a function  of  the  first  pass  rut  depth 

and  the  other  first  pass  parameters; 


M 

Cl 


instantaneous  first  pass  sinkage 
vertical  force  over  cone  index 


foot  print  area  over  cone  index 


I = foot  print  length 


The  variation  observed  in  all  tests  was  approximately  linear  and  hence  the 
rut  depth  of  an  nth  pass  is: 


Z = Z 
Rn  R1 


(n-1) 
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Incremental  Rut  Depth  - Direct  and  Indirect 


Data  were  collected  for  tests  condvicted  on  three  adjacent  paths. 

Rut  depths  were  examined  and  plots  made  to  determine  the  history  of  a 
particular  path  as  a function  of  direct,  load  applied  to  the  path,  and  indirect, 
load  applied  to  the  adjacent  path,  passes.  It  is  assumed  that  adjacent  lanes 
are  influenced  but  that  a separation  of  two  lanes,  two  tire  widths,  generates 
negligible  effect. 

Rut  depth  per  pass  for  the  tests  conducted  can  be  tabulated  as  shown 


below. 

R^ 

B^ 

B, 

P 

P/CI 

Test  Number 

D 

I 

D 

I 

A-0004-3 

0.  10 

-0.02 

0.  11 

-0.05 

750 

15 

A-0005-3 

0.25 

-0.00 

0.75 

-0.  35 

1060 

25 

A-0008-3 

0.06 

-0.02 

0.24 

-0.  12 

1150 

24 

A-0009-3 

0,  10 

-0.08 

0.28 

-0.02 

1160 

31 

The  terms  R^^  and  Rj 

are  t?ie 

rut  depth  added  per  pass  to  direct 

rolling  and  indirect  rolling.  The  terms  and  are  for  braked  response. 

The  data  were  collected  for  one  slip  value  and  there  was  no  apparent 
means  of  correlating  the  slip  to  the  incremental  rut  depth  values.  It 
appears  as  though  the  direct  braking  does  create  a rut  depth  change  pro- 
portional to  braked  sinkages,  about  3 to  1,  for  the  low  load  to  cone  index 
ratio.  However,  the  one  test  at  a low  load  for  the  same  slip  generates 
braked  rut  and  rolling  rut  of  the  same  magnitude.  Therefore,  until  tests 
are  conducted  at  a different  value  of  slip,  it  will  be  necessary  to  remember 
that  the  braked  incremental  rut  depth  is  directly  related  to  35  percent  slip. 

For  computational  purposes,  a graph  has  been  constructed  to  provide 
a guideline  for  analysis. 


o 

O Bj/B^ 
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Define  Multipass  History 

It  is  necessary  to  list  the  passes  in  order  as  they  occur.  The  history 
would  be  of  the  form; 


Pass  1 


IZ2£ 

Rolling 

Braked 

Braked 

Rolling 

Braked 

Braked 

Braked 


n Braked  ; 

The  rolling  path  effects  are  reasonably  well  defined,  but  it 
must  be  remembered  that  the  braked  effects  are  valid  only  for  35  percent 

slip.  It  is  assumed  that  the  paths  are  multiples  of  tire  width  so  that  path 
1 is  directly  down  the  middle  of  the  runway  and  all  others  are  offset  by  an 
integer  multiple  of  the  tire  width  distance. 


p , 

mi. 
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Calculate  Cumulative  Effects 

Cumulative  rut  depth  is  calculated  by  generating  a running  total  of 
the  effects  that  exist  on  each  path.  It  is  necessary  to  do  this  because 
although  over  100  passes  the  cumulative  results  may  be  acceptable,  the 
rut  depth  at  any  intermediate  pass  may  not  be.  This  is  true  because  indirect 
effects  apparently  raise  tlie  soil  and  reduce  the  rut.  Consequently,  many 
direct  passes  over  one  lane  could  generate  a rut  that  exceeded  the  criteria 
at  the  tenth  pass,  but  the  indirect  effects  of  the  11th  and  12th  pass  would 
make  a 12  pass  analysis  appear  acceptable  if  we  were  only  concerned  with 
the  cumulative  value  and  not  running  totals. 

The  final  step  is  to  tabulate  the  incremental  effects  with  each  pass. 

The  running  total  would  appear  as  shown  in  Figure  49. 


I 


G.  SUMMARY 

A procedure  has  been  described  to  calculate  the  effects  of  multipass 
on  soil  surfaces.  The  data  used  restricts  the  einalysis  to  the  following; 

1)  Cohesive  soil 

2)  Braked  drag  of  35  percent  slip 

3)  Rut  depths  of  approximately  three  inches. 

Additionally,  the  following  intuitive  restriction  is  necessary.  The  data 
indicated  that  indeed  the  soil  of  an  adjacent  lane  can  be  raised  above  its 
original  surface  by  a small  amount.  It  is  not  believed  that  this  can  exceed 
a very  small  eimount.  If  calculations  indicate  that  the  indirect  effects  want 
to  raise  the  soil  level  above  the  original  ground  level,  the  original  value 
should  be  used.  Cumulative  depression  is  reasonable,  cumulative  "heaving" 
above  ground  level  is  not. 
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I. 


Input  D^tu 

Tire  Type  7:00-6 

Dcncction  35% 

Contact  Area  43.4" 

Footprint  Length  9*0" 

First  Pass  Data 

Rolling  Sinkage  Z =0.85  Rolling  Rut  Depth  0.32 

Brak^::d  Sinkage  Z = 1.79  Braked  Rut  Depth  1.  16 

D 

Incremental  Rut  Depth  Data 

Incremental  Rolling  Rut  Depth  = 0.  30  = Rj^ 

Incremental  Braked  Rut  Depth  =0.90  = Bj^ 

Incremental  Rut  Ratios 

R B 

TT  = 0.25  =0.50 


Bj  = 0.  45 


Multipass  History  - Incremental  Ruts 


Pass  Zone 
No-  No, 


Zones 


1 


-0.08  0.32 

1.16  -0.58 

-0.08 

-0.  58 
-0.  08 
0.30  -0.08 


ulative  Ruts 


-0.08 

0.  32 

1.08 

-0.26 

1.08 

-0.26 

-0.26 

1.00 

-0.26 

-0. 26 

1.00 

-0.84 

-0.26 

1.00 

-0.92 

-0.34 

1.30 

-1.00 

-0.58 

-0.08 


Comments; 
First  Pass  on  1 
First  Pass  on  2 
First  Pass  on  5 
First  Pass  on  4 
First  Pass  on  3 


).08 

>.08 

wBm 

).  16 

).  16 

1.00 

-0.26 

1.30 

-0.  34 

1.30 

-0.  34 

M- 


Figure  49.  Sample  Cumulative  Rut  Calculations, 


SECTION  IV 

STATIC  STARTUP  FORCE  ANALYSIS 

A.  INTRODUCTION 

Another  aspect  of  tire /soil  interaction  is  that  of  the  difference  between 
start-up  forces  and  rolling  drag  forces  for  a given  tire  in  a particular  soil. 

It  has  been  recognized  that  the  difference  does  exist,  but  no  quantitative  means 
of  predicting  this  is  available.  Consequently,  existing  data  were  examined  to 
determine  how  well  it  could  be  predicted  using  rolling  tire  prediction 
parameters. 


B.  DATA  AVAILABLE 

The  sources  for  start-up  force  data  are  listed  in  Table  I on  the 
following  page.  These  were  selected  because  they  were  the  only  single  wheel 
tests  conducted  where  initial  drawbar  pull  data  existed  along  with  peak 
rolling  forces  cind  average  rolling  force.  All  tests  were  conducted  in  a hea\’y 
clay  soil  (buckshot). 


I\V 


1 
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C.  DATA  PRESENTATION 

The  primary  concern  of  previous  efforts  has  been  the  development  of 

rolling  or  braked  drag  ratios  that  could  be  predicted  based  upon  tire/soil 

parameter  combinations.  Similarly,  then,  static  start-up  drag  ratios  should 

also  be  calculated  for  the  same  tire/soil  parameter  combinations  if  possible, 

R 

Consequently,  the  static  start-up  drag  rario  is  defined  as  _S  , where  R 

P ^ 

is  the  start-up  drag  force  and  P the  applied  single  wheel  vertical  load.  It 
was  assumed  that  just  as  rolling  drag  was  dependent  upon  the  ratio  of  applied 
stress  to  cone  index  (a /Cl),  so  should  the  static  start-up  drag. 

Figure  50  presents  all  of  the  data  available  for  the  Z5.00  x 28,  30  ply 
tire.  Two  curves  are  shown  for  both  free  rolling  drag  ratio  and  start-up 
drag  ratio  since  both  first  pass  and  twentieth  pass  data  were  available.  The 
start-up  ratios  show  considerable  scatter  and  the  free  rolling  valuer  appeared 
to  be  low.  Magnitudes  of  the  free  rolling  drag  ratio  were  computed  using  the 


TABLE  I 
DATA  SOURCES 


Reference  Repor 


AFFDL-TR-66-43 
Part  m 

AFFDL-TR-66-43 
Part  Vn 


AFFDL-TR-66-43 
Part  Vm 


AFFDL-TR-66-43 
Part  IX 


AFFDL-TR-66-43 
Part  X 


Tire  Type 

Inflation 

Pressure 

56x16,  24  Ply 

100  psi 

25x28,  30  Ply 

50  psi 

25x28,  30  Ply 

100  psi 

56x16,  32  Ply 

100  psi 

25x28,  30  Ply 

100  psi 

17x16,  12  Ply 

100  psi 

34x9.9, 14  Ply 

100  psi 

25x28,  30  Ply 

25  psi 

25  X 28,  30  Ply 

40  psi 

25x28,  30  Ply 

60  psi 

25x28,  30  Ply 

80  psi 

*25x28,  30  Ply 

100  psi 

25x28,  30  Ply 

50  psi 

12.0 

9.2 


Applied 
Vertical  Load 


35,000 


35.000 

60.000 

25,000 

25,000 

25,000 

25,000 

25,000 

25,000 

25,000 

25,000 

25.000 

35.000 


* This  test  was  not  used  since  lane  trafficked  had  been  used  before 
although  "essentially  undamaged". 


Drag  Ratio  (R/P) 


0.20 


Start-up  First  Pass 


0.18  - 

0.16  ■ 
0. 14  - 
0. 12  - 
0.10  - 
0.08  ■ 
0.06  ■ 
0.04  - 
0.02  - 

0.0  + 

0 


, j 1 "I  "I 

0.  1 0.2  0.3  0.4 


g 

Cl 

LOAD  STRENGTH  RATIO 


Figure  50.  Drag  Ratio  Comparisons  versus  Load  Strength  Ratio, 
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available  predictive  equations  and  plotted  as  shown  on  Figure  5j  . The 
load  strength  ratios  (a  /Cl)  were  calculated  using  the  contact  areas  given 
in  the  referenced  reports  for  the  particular  tire  at  the  given  tire  pressure. 
Since  there  was  such  a significant  difference  between  measured  and  calculated, 
the  peak  rolling  forces  were  plotted  as  also  shown.  The  results  of  examining 
the  "peak"  rolling  and  "average"  rolling  ratios  is  that  it  is  apparent  that  the 
current  predictive  equations  reflect  the  average  force  ratios  for  the  two. 

The  start-up  force  ratio  can  be  approximated  by  using  the  slope  of  the 
cadculated  rolling  drag  ratio  and  adding  ein  intercept  magnitude  of  0.  1. 
Therefore: 


R 

P 


+ 0. 1 


is  the  gross  approximation  of  the  data. 

It  is  possible  to  conduct  further  analysis  on  the  data  but  it  is  not 
believed  to  be  justified.  The  data  points  seem  to  be  functions  of  cone  index, 
applied  load  and  tire  pressure  although  plotted  against  a/CI  which  contains 
the  same  parameters.  K only  a given  applied  load  and  cone  index  is 
examined  there  are  four  data  points  that  can  be  used.  These  indicate. 

Figure  52  , that  there  may  be  a tendency  to  follow  the  linear  approximation, 
or  that  there  may  be  a curve  in  the  response  observed. 

The  linear  plot  could  be  an  approximation  because  of  the  data  scatter. 
The  curve  can  also  be  justied.  At  the  lower  a/CI  values  the  tire  is  soft 
relative  to  the  soil.  If  the  tire  becomes  softer,  it  ultimately  approaches  a 
flat  tire  which  indeed  has  a large  start-up  force.  At  the  higher  a/CI  value 
the  tire  has  greater  stiffness  than  the  soil  and  in  the  lim.t  acts  as  a thin 
sharp  tire  sinking  deeply  into  a soft  soil.  Both  extremes  could  suggest  that 
there  is  some  optimum  value  of  applied  force,  tire  pressure  and  soil  index 
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Drag  Ratio  (R/P) 


to  minimize  atart-up  drag.  However,  the  data  are  too  limited  to  find  the 
true  nature  of  the  atart-up  approximation. 


r 


i 


Other  attempta  were  made  to  iaolate  the  effecta  of  tire  presaure  for 
a given  cone  index,  applied  force  for  a given  preaaure,  and  other  combinations 
of  parametera.  However,  the  data  are  too  limited  for  any  rigorous  analysis. 

D.  SUMMARY 

The  data  examined  indicate  that  the  start-up  drag  force  is  greater  than 
the  measured  or  predicted  rolling  drag  force  for  any  given  load  strength 
ratio.  The  number  of  data  points  available  and  the  number  of  parameters  to 
be  isolated  are  not  compatible  in  terms  of  finding  exact  relations.  The  trend 
of  start-up  drag  can  be  seen,  but  how  it  can  be  more  exactly  predicted  cannot 
be  determined.  For  preliminary  criteria,  if  necessary,  the  start-up  drag  ratio 
is: 


R 

P 


+ 0.1 


recognizing  that  this  was  found  from  tests  conducted  in  buckshot  clay  with 
the  pressures  of  Z5  to  100  psi,  vertical  forces  of  ?.5,  000  to  60,  000  pounds, 
and  soil  strengths  of  CBR  3.9  to  CBR  12.  0, 


SECTION  V 

TIRE /SOIL  ROUGHNESS  INTERACTION  STUDIES 


A.  INTRODUCTION 

Another  area  of  the  research  effort  to  be  considered  was  tlie  investi- 
gation of  surface  roughness  effects  upon  tire  and  landing  gear  loads.  There 
have  been  many  research  efforts  directed  toward  solving  the  problem  of  tire 
and  soil  interaction  with  roughness  present.  Unfortunately,  there  have  only 
been  a few  analyses  where  flexible  tires,  soft  soil  and  roughness  have  all 
been  considered  simultaneously.  The  purpose  of  the  analytical  representation 
selected  to  compute  such  an  interaction,  was  that  it  should  reflect  tire 
deflection,  soil  compression,  drawbar  pull  or  driving  torque,  surface  slip 
and  the  nature  of  the  soil,  whether  cohesive  or  cohesionless. 

There  have  been  other  representations  similar  to  that  described,  but 
it  is  believed  that  the  current  model  has  greater  flexibility  and  potential 
capability  to  predict  the  desired  response. 

B.  BACKGROUND 

The  tire  as  it  rolls  through  the  soil  deforms  and  depresses  the  soil. 

At  the  interface  between  the  two,  shear  and  normal  forces  act  to  create  an 
equilibrium  between  horizontal,  vertical  and  angvilar  forces  and  inertial 
accelerations.  The  problem  classically  has  been  to  evolve  a set  of  equations 
that  predicts  evaluation  of  the  forces  at  the  interface  which  will  generate 
the  measured  forces  at  the  wheel  hub  and  the  soil  instantaneous  sinkage  euid 
subsequent  rut  depth.  The  problem  is  further  complicated  if  we  permit  the 
tire  to  roll  over  large  obstacles  which  cause  localized  forces  and  deformations 
significantly  different  from  what  would  be  called  "steady- state  rolling". 

The  equations  that  have  been  used  to  present  the  tire/ soil  response 
can  be  best  explained  by  using  ein  auialytical  model.  The  model  consists  of 
several  elements  selected  to  best  duplicate  the  tire  amd  soil  without  being 
excessively  complicated. 
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The  tire  is  represented  by  a series  of  rotating  spring  mass-damper 
systems.  Figure  53  shows  a tire  segment  which  represents  a nine-degree 
portion  of  the  tire.  The  wheel  hub  is  the  attachment  point  for  all  radial 
segments,  and  is  the  location  of  the  application  of  all  vehicle  inputs.  That 
is,  at  the  hub  w e have  the  vehicle  initial  vertical  and  horizontal  motions 
and  forces.  The  tire  "tread"  is  represented  by  a mass  which  is  separated 
from  the  hub  by  an  elastic  and  viscous  element.  Usually  these  are  linearly 
elastic  elements,  but  for  the  computer  routine  evolved,  this  is  not  mandatory. 
The  values  of  the  stiffness,  damping  and  mass  are  inputs  to  the  program. 

At  the  tire  inertial  element  is  shown  the  soil  element.  The  auialytical 
representation  was  selected  because  it  was  desirable  to  have  damped 
inertial  response,  the  top  elements,  and  some  representation  of  permanent 
deformation,  the  lower  element.  It  is  known  that  the  drag  force  of  the  tire 
has  a peak  at  some  speed  and  then  drops  off.  The  model  used  permits  this 
to  occur  while  allowing  a calculation  of  both  instantaneous  sinkage  and  x*ut 
depth. 

As  shown,  the  figure  indicates  one  tire  segment  phishing  down  upon 
one  soil  element.  In  order  to  have  a more  valid  representation  of  the  tire, 
a series  of  these  segments  are  used.  Figure  54  indicates  11  elements  used 
to  represent  the  lower  portion  of  the  tire.  These  rotate  at  an  angular  velocity 
and  translate  forward  as  dictated  by  the  velocities  of  the  hub.  As  the  wheel 
"rolls",  each  tire  segment  contacts  a soil  element  and  compresses  it.  If 
at  the  surface  we  have  a normal  and  tangential  force,  it  is  possible  to  write 
the  equation  of  equilibrium  as  shown  in  Figure  55. 

The  free  body  diagram  indicates  that  a moment  is  applied  at  the  hub. 

This  can  be  either  positive  or  negative  depending  upon  whether  or  not  the 
tire  is  torque  driven  or  braked.  Each  tire  segment  contributes  to  the 
summation  of  forces  and  moments  as  it  rotates  from  a starting  position  through 
an  arc  of  90  degrees. 
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Figure  53.  Tire/Soil/Roughness  Model. 
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The  soil  surface  is  represented  by  a series  of  line  segments  defined 
by  vertical  position  as  functions  of  displacements  along  the  "runway”.  The 
points  are  input  into  the  program  as  ordered  pairs. 

The  calculation  of  resulting  tire/soil  forces  is  made  by  following  the 
position  of  the  tire  "tread"  and  determining  where  and  when  the  intersection 
of  tire  amd  soil  contour  begins.  When  this  occurs,  the  integration  of  the 
differential  equations  begins  and  the  motion  of  the  hub  results.  The  outputs 
from  the  program  are  hub  motions  and  forces,  tire /soil  interface  forces, 
instantcineous  sinkage  and  rut  depth. 

The  Soil  Element 

The  soil  is  modeled  by  having  a spring -mass -damper  system  resting 
on  top  of  a viscous  element.  This  configuration  was  chosen  because  of  its 
capability  to  produce  a "resonant"  force.  Additionally,  the  two  dependent 
parameter  model  periruts  instantaneous  sinkage  as  well  as  rut  depth 
prediction. 

The  model  is  shown  below  with  a positive  applied  force. 
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The  applied  force,  F,  is  developed  at  the  tire  tread  inter- 
face. The  upper  stiffness  and  damping  elements  are  K_  and  C , respectively. 

o S 

The  soil  mass,  that  beneath  the  tire  tread,  is  M , and  the  "subsoil”  damping 

is  characterized  by  C . The  soil  motion  is  all  referenced  to  an  inertial 

00 

ground. 

At  the  applied  force, 

E F = Ma  = MZ^ 
since  no  mass  is  present 
E F = 0 

• • 

F,  + - z^)  + - Z)  =0 

1 S S T S S T 

or, 

= ‘'s  <^T  - ^S>  + S '^T  - ^S>  • 


The  force  supplied  by  the  tire  equals  the  resistance  developed  by  the  "upper" 
soil. 


or. 


For  the  "soil"  element. 


=K^  (Z^  - ZJ  + (Z^  - Z^)  - C^^Z^ 

S S S T S S ' T s'  SS  S 
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These  are  the  equations  necessary  to  define  the  response  of  the  soil.  The 
only  remaining  aspect  is  that  a logiceil  control  must  be  added  such  that  the 
force  at  the  tire  tread  cannot  be  a tensile  force.  The  soil  cannot  pull  the 
tread.  When  the  force  gets  to  zero,  all  soil  forced  response  must  stop  and 
transient  response  follow. 


The  soil  response  starts  whenever  the  tread  contacts  the  soil. 

At  small  angles  of  tire  impact,  when  6 is  small,  the  model  is  intuitively 
unrealistic  since  it  is  likely  that  a cohesionless  soil  would  be  pushed  out  of 
the  way  rather  than  compressed.  However,  in  order  for  the  equations  of  soil 
and  tread  to  be  consistent,  the  same  action/ reaction  force  and  displacements 
for  both,  it  is  necessary  to  model  the  interaction  as  though  only  compression 
occurs. 

Later  it  will  be  shown  that  the  normal  force  component  that  acts  on 
the  tire  will  be  used  to  calciilate  the  frictional  forces  developed. 

Tire  Tread 

A portion  of  the  tire  is  represented  by  the  tire  "tread”  at  the  interface 
between  tire  and  soil. 


The  tread  is  acted  upon  by  the  soil,  a normal  and  frictional  force,  the 
inertial  effects  of  the  rotation,  and  the  elasticity  and  damping  of  the  tire. 
The  forces  are  shown  with  the  tire  accelerating  positively.  Summing 
forces  in  the  vertical  direction. 


It  is  assumed  that  the  traction  across  the  tread  is  transferred  directly  without 


any  moment.  Therefore  T = T and  the  equation  becomes 

X M 


. 2 


T ~ ^ sin9  + M^R  0 cos6  - R0  sin  9 - sin  9 - 


The  forces  due  to  the  angular  motion  of  the  wheel  can  also  he 
eliminated  since  they  are  internally  accounted  for.  That  is,  the  centrifugal 
force  of  the  element  is  balanced  by  a component  of  the  tire  reaction  F^. 

A similar  component  exists  on  the  other  side  of  the  wheel  hub,  and  hence 
the  net  effect  upon  wheel  response  is  negligible.  The  component  due  to 
angular  acceleration  is  accounted  for  in  that  no  net  force  on  hub  is  created 
and  the  moment  generated  is  contained  within  the  product  of  the  mass 
moment  of  inertia  of  the  entire  wheel  multiplied  by  the  wheel  angular 
acceleration.  Therefore,  the  tread  equation  reduces  to 


M„Z  ^ = N sin  9 - F^  sin  9 - 
T T T T 


The  normal  force,  N,  is  the  normal  of  the  soil  force  F_  , 

DX 


hence 


N = F sin  9 in  magnitude. 
1 SI 


However,  since  F^^  was  derived  for  a positive  displacement  of  the  soil  and 


the  force  acting  upon  the  tire  is  equal  and  opposite  to  that  acting  on  the  soil. 


N = -F  sin  9 , is  the  correct  equation  for  the  free  body  of  the  tread. 
X ox 


At  the  surface  of  the  tire  a frictional  force  is  developed  which 
is  carried  into  the  elastic  side  wall.  The  equation  used  for  this  was; 


i#" 


T = A (o  tan  <p  - C)  \i  (s) 

Ca 

where  T is  the  frictional  forces  in  pounds 
A is  the  tread  area,  square  inches 
a is  the  normal  stress,  pounds/square  inch 
tain  0 is  the  angle  of  internal  soil  friction,  radians 
C is  the  soil  ''ohesion,  pounds /square  inch. 

/I  (s)  is  the  cube  root  of  the  slip. 

The  equation  was  developed  in  Reference  3 and  provides  a means  of  cal- 
culating the  surface  force  as  a function  of  both  soil  properties,  and  local 
slip. 

The  area  of  the  tread  is  calculated  assuming  that  each  tread 
segment  represents  nine  degrees  of  arc.  This,  with  the  tire  diameter  and 
section  width,  is  sufficient  to  calculate  a tread  area.  Specifically,  for 
nine  degrees, 

=B^  (.  15708)  RW 

where  is  the  section  width,  and  RW  is  the  tire  radius. 


Slip  is  defined  as 


SLIP  = 


V - V 
w a 


where  V is  the  peripheral  speed  of  the  wheel,  and  V is  the  horizontal 
w ■ a 

speed.  However,  for  the  portion  of  the  tire  which  is  touching  a slope, 
the  slip  must  be  redefined  to  account  for  the  components  along  the  slope. 


Slip  is  now  the  difference  between  peripheral  speed  and  the 
component  of  the  horizontal  speed  divided  by  the  peripheral  speed 

SLIP  = R9.  :.X,  sin, 6 

R9 

This  reverts  to  the  original  definition  for  tangency  at  the  bottom  of  the 
wheel. 

If  the  peripheral  speed  is  greater  than  the  translational  com- 
ponent, the  slip  is  positive  as  a powered  wheel,  whereas  negative  slip  infers 
braking. 

The  frictional  force,  T , originally  defined  in  the  tire  tread 
free  body  is  shown  as  a positive  force  upward.  This  is  consistent  with  a 
positive  slip.  If  the  peripheral  velocity  exceeds  the  translational  component, 
slip  is  positive,  the  wheel  slips  over  the  soil  and  a friction  force  is  developed 
that  acts  to  retard  the  motion  in  a positive  direction. 

Tire  Elastic -Viscous  Element 

The  tire  "side"  wall  is  an  elastic  and  viscous  element  in  parallel. 
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It  is  assumed  that  the  friction  force,  T , can  be  carried  by  the  bending 
stiffness  of  the  side  wall  element  without  deformation.  Therefore,  T is 
balanced  by  T and  at  the  top  of  the  element,  the  wheel  hub.  The  element 
is  an  axial  force  element  and  the  force  is: 


Ft  = (6^  - 6^)  + (6^  - 6^) 


where  is  the  stiffness  in  pounds /inch 

is  the  damping  in  pounds/inch/second 

is  the  axial  deformation  at  the  tire  tread 

is  the  axial  deformation  of  the  hub. 

w 


Since 


= 6 ^ sin  e , and 


^'W  = ® 


Wheel  Hub 


. (K^  (Z^  - Z^)  + C^  (Z^  - Z^)  ) /sin  6 . 


The  free  body  diagram  for  the  wheel  hub  is  shown  below. 


where  DP  is  the  drawbar  pull, 

W is  the  weight  of  the  payload 

w 

Wp  is  the  weight  of  the  payload 
M is  an  applied  torque  (positive  or  negative) 

P is  an  applied  vertical  force. 

The  equations  of  motion  are: 

« • 

(Mp  + M^)  Z ^ = T cos  9 + sin  9 - - P - W 

• • 

(Mp  + M^)  = T sin  9 - cos  9 + DP 

V 9 ~ ~ ’ 


The  only  variable  not  previously  defined  is  M^.  This  is  the  moment  gen- 


erated by  the  friction  it  the  soil  and  carried  through  the  side  wall.  It  is 
therefore, 


= T . r 
T 


where 


r = 


“ ^ " ^T  ’ 


the  undeformed  radius  RW,  plus  the  extension  of  the  hub,  minus  compression 
at  the  tread.  Since 


6 ^ = Z^/sin  e , 6 ^ = Z^/  sin  0 , 


then 


r = RW  + (Z,„  - Z)  / sin  0 

W 1 


Logical  Controls 


The  equations  developed  are  necessary  to  calculate  the  force  and 
moment  response  of  a soft  tire  as  it  compresses  a soft  soil.  However,  the 
equations  reflect  the  response  as  thou^  at  some  initial  time  tire  and 
soil  are  touching  with  known  initial  conditions.  This  is  not  true  in  that  we 
wish  to  have  a tire  element  rotate  until  it  touches  the  soil,  and  then  the 
response  is  calculated  for  whatever  conditions  existed  at  that  moment. 


Recognition  of  tire-soil  impact  is  accomplished  by  •'omputing  the 
difference  in  elevation  between  the  tread  and  soil  at  any  particular  time. 


-*-X 
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c. 


ROLLING  TIRE  COMPUTER  PROGRAM  AND  SAMPLE  RESULTS 


The  equations  evolved  were  programmed  using  the  MIMIC  computer 
program  detailed  in  Reference  19.  A complete  listing  of  the  program  is 
shown  in  Appendix  B.  Many  comment  cards  have  been  included  in  the 
printout  to  indicate  the  appropriate  steps  being  taken  for  each  block  of 
cards. 

A sample  run  was  made  using  the  following  input  conditions. 


Initial  horizontal  velocity 

240. 

inches  per  second 

Initial  vertical  velocity 

0. 

inches  per  second 

Initial  vertical  applied  force 

0. 

pounds 

Drawbar  pull 

500. 

pounds 

Braking  moment 

0. 

inch  pounds 

Payload  weight 

500. 

pounds 

Wheel  weight 

500. 

pounds 

Wheel  height  above  ground  reference  30. 

inches 

Tire  radius 

25. 

inches 

Wheel  moment  of  inertia 

500. 

pound-inch  -seconds  square 

Tire  segment  weight 

5. 

pounds 

Tire  Segment  stiffness 

500. 

pounds 

Tire  segment  damping 

5. 

pounds  per  inch  per  second 

Tire  width 

20. 

inches 

Weight  of  soil  segment 

20. 

poTUids 

Soil  segment  stiffness 

2000. 

pounds  per  inch 

Soil  segment  damping 

0. 

pounds  per  inch  per  second 

Subsoil  daunping 

20. 

pounds  per  inch  per  second 

Height  of  soil  above  reference 

6.5 

inches 

Soil  cohesion 

1. 

pound  per  square  inch 

Internal  friction  angle 


0. 016  radian 


The  numbers  used  are  generally  indicative  of  realistic  values  although 
some  of  them  are  strictly  estimates.  The  soil  model  characteristics  are 
not  knowns  since  we  do  not  know  how  a soil  "segment"  reacts.  This  doesn't, 
however,  eliminate  the  use  of  real  world  value  for  cohesion  and  integral 
friction  at  the  tire  tread/soil  interface.  Consequently,  those  values  are  for 
a low  strength  clay  soil. 

The  computed  drag  and  slip  are  shown  on  Figure  56.  The  pilots 
indicate  that  in  free  falling  for  one-tenth  of  a second,  transient  response  is 
still  present  alter  0.  65  seconds  of  motion.  Assuming  that  steady  state  drag 
will  be  about  100  pounds,  the  drag  ratio  would  be  0.1.  The  sinkage  computed 
for  the  same  impact  is  approximately  2 inches  for  a 50 -inch  diameter  wheel. 
Therefore,  Z/D  = 0.  04  and  this  corresponds  to  a rolling  drag  ratio  of  0. 17 
for  CBR  6 soil,  v_-onsequently  the  drag  is  at  least  the  right  order  of 
magnitude  for  a rather  hypothetical  wheel. 

The  computation  shown  required  approximately  five  minutes  of  com- 
puter time  using  a CDC  6600  series  computer.  If  the  routine  had  been  inm 
long  enough  for  a steady  state  solution  and  then  a "roughness"  introduced, 
it  is  apparent  that  significant  computer  time  would  be  required.  Additional 
study  is  required  to  more  thoroughly  evaluate  the  model  and  establish 
whether  or  not  sufficiently  accurate  results  can  be  established  with  a single 
repre  sentatiLon. 

D.  ROUGHNESS  INTERACTION  STUDIES  SUMMARY 

A computer  program  has  been  developed  which  is  designed  to  calcvilate 
the  effects  of  soil  surface  roughness  upon  the  sinkage  and  drag  of  a pneumatic 
tire.  The  program  has  the  capability  to  utilize  true  soil  properties  and 
realistic  tire  and  wheel  characteristics  and  applied  loads,  and  surface 
roughness.  The  one  aspect  remaining  to  be  solved  is  that  of  defining  tlie 
analytical  characteristics  of  soil  euid  tire  stiffness  eind  damping.  The 
program  was  developed  in  order  to  use  as  input  existing  test  data  for  rolling 


tires.  Since  the  test  data  would  reflect  the  input  and  outputs  for  given  soils 
and  tires,  the  model  coefficients  would  be  adjusted  to  match  the  measured 
output  of  drag,  instantaneous  sinkage  and  rut  depth.  The  program  was 
developed  to  the  point  of  generating  the  desired  outputs,  but  it  was  not 
possible  to  conclude  the  effort  with  test  data  which  are  available.  Hopefully, 
in  the  future  such  steps  will  be  taken. 


SECTION  VI 

ADDITIONAL  EFFORT  IN  TIRE/SOIL  INTERACTION 


A.  TRAINING  SESSION  - SHORT  COURSE  ON  LANDING  GEAR /SOIL 

INTERACTION  AND  FLOTATION  SYSTEM  DESIGN 

A short  course  on  landing  gear/soil  interaction  and  flotation  design 
was  held  at  the  University  of  Dayton  on  September  18  and  19,  1974.  The 
course  was  designed  to  disseminate  the  latest  information  available  on 
landing  gear/soil  interaction  phenomena  and  associated  aircraft  design 
Bind  operation  criteria  for  soil  runways  to  government  and  industry 
personnel,  A total  of  21  persons,  representing  airframe  manufacturers, 
aircraft  tire  manufacturers,  governmental  agencies,  and  one  university 
were  present. 

Three  volumes  of  material  were  provided  to  all  participants  along 
with  video  tape  explanations  and  workshop  problem  sessions.  The  areas 
discussed  were; 

- Aircraft  tire/soil  sinkage  and  drag  performance 

- Multiwheel  sinkage  and  drag  performance 

- High  speed  sinkage  and  drag  performance 

- Braked  tire  sinkage  and  drag  performance 

- Turning  performance  on  soils 

- Takeoff  and  landing  performance  on  soils 

A complete  set  of  instructional  materials  can  be  obtained  for  the  cost 
of  reproduction  upon  request  to: 

Dr,  David  C.  Kraft 
Dean  of  Engineering 
University  of  Dayton 
Dayton,  Ohio  45469 
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SECTION  vn 

DEVELOPMENT  OF  A FULL-SCALE  TEST  PROGRAM 


A.  INTRODUCTION 

The  following  proposed  test  plan  outlines  a series  of  tests  to  be  con- 
ducted using  a full-scale  aircraft  operating  on  a soil  runway.  The  type  of 
aircraft  is  not  specified,  but  the  parameters  to  be  controlled  and  measured 
should  be  applicable  to  any  aircraft  having  a medivim  transport  classification, 

B.  PURPOSE 

Significant  advances  have  been  made  in  the  last  several  years  in  the 
area  of  predicting  tire/ soil  response  of  aircraft  tires  while  free  rolling, 
turning,  or  being  braked.  The  majority  of  the  data  has  been  the  results  of 
tests  conducted  in  laboratories  using  loads  and  tires  that  are  not  those  seen 
in  operational  conditions.  Although  the  prediction  equations  are  based  upon 
particular  tire /soil  parameters  which  correlate  with  a great  amount  of  data 
collected  for  many  tire  types  eind  sizes  and  several  soil  types  and  strengths; 
still  it  remains  to  conduct  full-scale  tests. 

The  tests  proposed  are  concerned  with  the  taixi  and  turning  aspects  of 
ground  maneuvering.  Turning  is  a function  of  free  rolling  drag,  braked  drag 
and  side  forces  developed  on  both  nose  wheel  and  main  gear.  Since  the  side 
forces  are  functions  of  the  zero  steering  angle  forces,  it  is  necessary  to 
measure  "taxi"  loads  first,  before  proceeding  on  to  the  turning  mameuver. 

The  accomplishment  of  the  proposed  tests  would  accomplish  several 
objectives; 

1)  Permit  comparison  of  measured  forces  on  full-scale  wheel 
assemblies  under  operational  loads,  with  predicted  values  for 
free  rolling  and  braked  wheels. 


W- 


V 
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I 2)  Permit  comparisons  of  measured  and  calculated  rut  depths  for 

f 

I free  rolling  and  braked,  zero  steering  angle,  wheel  assemblies. 

I 

I 3)  Provide  guidance  in  updating  the  existing  routines  so  that  additional 

parameters  that  are  significant  can  be  isolated  and  incorporated. 


C.  TEST  PROGRAM 

There  are  many  ways  of  establishing  the  test  matrix  for  a particular 
test  program.  The  most  direct  approach  for  this  program  is  to  consider  the 
various  means  of  performing  over  the  ground  maneuvers.  Specifically,  there 
are  three  ways  to  maneuver  over  the  ground  at  low  airspeeds.  These  are; 

1)  differential  engine  thrust,  2)  differential  braking,  and  3)  nose  wheel 
steering.  These  can  be  tabulated  as  shown  below. 


Brakes 

Thrust 

(Symmetrical) 

Thrust  (Unsymmetrical)  | 

Trail 

Steering 

Trail 

Steering 

None 

1 

3 

2 

4 

Symmetrical 

5 

7 

6 

8 

Differential 

9 

11 

10 

12 

This  suggests  the  hierarchy  of  testing  required  to  accurately  isolate 
the  many  variables  present.  For  example,  block  number  1;  syinmetrical 
thrust,  no  steering  and  no  brake,  is  a description  of  a free -rolling  taxi 
maneuver.  If  a vehicle  were  tcixiiing  at  constant  thrust  and  velocity  under 
these  conditions,  any  measurements  taken  on  the  landing  gear  mounts  would 
be  indicative  of  free  rolling  drag  as  modified  by  any  multiwheel  effects. 
Therefore,  for  the  particular  wheel  and  tire  configuration,  load  distribution 
and  soil,  it  is  possible  to  compare  the  forces  developed  with  predicted 
values  from  existing  equations. 
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The  next  added  complexity  could  then  be  that  of  braking  as  listed  in 
block  number  5.  Symmetrical  braking  of  the  main  gear  does  not  cause  any 
turning  (theoretically)  but  does  cause  increased  drag  loads  as  a function  of 
wheel  slip.  No  turn  angle  or  steering  angles  would  be  introduced  and  it  would 
be  possible  to  collect  data  for  braked  wheel  response. 

All  other  cells  of  the  table  do  introduce  some  turning  and  hence 
generate  a turn  angle  for  all  wheel  assemblies.  The  question  is  how  to  best 
conduct  the  testing  to  minimize  any  multiple  parameter  effects  and  maximize 
the  effects  of  particular  parameters.  The  next  test  block  would  therefore  be 
number  3 which  introduces  nosewheel  steering  as  the  means  of  turning. 

By  using  only  nose  wheel  steering  with  free  rolling  wheels  and  constant 
thrust,  the  vehicle  has  to  reach  an  equilibrium  condition  where  there  will 
be  a specific  turn  angle  for  each  wheel  assembly  as  dictated  by  the  free 
rolling  side  and  drag  forces.  Additionally, at  that  condition  a turned  sinkage 
and  rut  depth  will  be  generated.  The  effects  of  thrust  and  braking  have  been 
eliminated. 

The  next  higher  order  of  complexity  would  be  possible  testing  with 
symmetrical  braking  as  identified  by  cell  number  7.  By  adding  symmetrical 
braking  it  is  possible  to  establish  the  effects  of  braked  slip  at  a turn  angle. 
The  effects  of  combined  braking  and  turning  on  tire  sinkage  and  drag  are  not 
known.  Since  the  previous  test  series  was  to  be  conducted  at  a turned  con- 
dition with  no  braking,  it  should  be  possible  to  estimate  how  the  sinkage  and 
drag  are  related  with  braking. 

This  completes  the  testing  required  in  that  any  other  combination  of 
thrust,  braking  or  steering  does  not  introduce  a new  environment  for  the 
wheel.  It  can  either  be  straight  ahead  rolling,  turned  rolling,  or  braked 
with  and  without  a turn  angle.  The  four  test  series  outlined  provide  data 
for  all  of  these  combinations.  Any  additional  testing  would  certainly  be 
beneficial  in  better  defining  the  effects  of  thrust  and  braking  upon  overall 
vehicle  response,  but  would  not,  theoretically,  alter  the  interaction  effects 
already  evaluated. 
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Further  refiiiement  of  the  test  program  is  possible  when  we  consider 
the  number  of  significant  parameters  and  the  possible  variations.  These 
are: 


Parameter 

Variations 

Number 

Gross  Weight 

Design  gross  weight  and  max  gross 
weight. 

2 

Taxi  Speed 

Minimum  and  maximum  of  Region  II 
(5  to  20  kts) 

2 

Tire  Deflection 

Operational  value  only  (32%) 

1 

Soil  Type 

Clay  or  Sand 

2 

Thrust 

Symmetrical  as  required  for 
constant  speed. 

1 

Steering 

Trail  position,  20,  40  and  60°. 

4 

Brakes 

Symmetrical,  15,  30,  and  45 
percent  slip. 

4 

The  steering  values  of  20,  40  and  60*^  were  selected  based  upon  known 
criteria  for  C-130E  turning.  The  braked  slip  values  are  established  to 
investigate  minimum  braking,  an  intermediate  value  which  theoretically 
generated  a braked  drag  ratio  of  50%  and  a maximum  value  which  should 
approach  the  0.  8 drag  ratio  required  by  specification. 

Using  the  above  values  generates  the  following  nvunber  of  tests  . 


'i 


Test  Type 

Control  Variables 

1 

5 ; 3 

7 

E 

Gross  Weight 

2* 

1 

2 2 

2 

Taxi  Speed 

2 

2 i 2 

2 

Soil  Type 

2 

2 2 

2 

inrust 

* 

Steering 

1 

1 3 

3 

Brakes 

1 

3 1 

3 

Product 

8 

24  ’ 24 

i 

72 

128 

* Number  of  variations 


1- 


ffe. 


Therefore,  1Z8  separate  passes  would  be  required.  Realistically  this  ixr, plies 
that  at  one  teat  site  64  taxi  tests  would  be  required  in  order  to  cover  a spectrum 
of  operational  conditions. 

Test  Matrix 

The  test  matrix  for  one  testing  site  would  be  as  shown. 

TABLE  n 

TEST  SITE  "A"  (SOIL  I) 


Test 

Number 

Gross 

Weight 

Velocity 

(Kts) 

Steering 

(Degrees) 

Braked 
(%  slip) 

Objectives 

101 

^1 

5 

Trail 

None 

Verify 

102 

^2 

5 

Trail 

None 

Free  Rolling 

103 

^1 

20 

Trail 

None 

Drag  and  Sinkage 

104 

^2 

20 

Trail 

None 

Prediction 

Techniques  - With 
no  Turning 

501 

Wi 

5 

Trail 

15 

Verify 

502 

^1 

5 

Trail 

30 

Braked  Rolling 

503 

^1 

5 

Trail 

45 

Response 

504 

^1 

20 

Trail 

15 

Prediction 

505 

^1 

20 

Trail 

30 

Techniques 

506 

^1 

20 

Trail 

45 

With  no 

507 

^2 

5 

Trail 

15 

Turning 

508 

^2 

5 

Trail 

30 
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TABLE  II  (Continued) 


Test 

Number 

Gross 

Weight 

(lbs) 

Velocity 

(Kts) 

Steering 

(Degrees) 

Braked 
(%  Slip) 

Objectives 

509 

^2 

5 

Trail 

45 

510 

^2 

20 

Trail 

15 

511 

^2 

20 

Trail 

30 

512 

^2 

20 

Trail 

45 

301 

^1 

5 

20 

None 

Establish 

302 

'"i 

5 

40 

None 

Effects  of 

303 

^1 

5 

60 

None 

Turn  Angle 

304 

^1 

20 

20 

None 

on 

305 

*1 

20 

40 

None 

Free  Rolling 

306 

^1 

20 

60 

None 

Wheels 

307 

^2 

5 

20 

None 

308 

^2 

5 

40 

None 

309 

^2 

5 

60 

None 

310 

^"2 

20 

20 

None 

311 

^2 

20 

40 

None 

312 

^2 

20 

60 

None 

701 

^1 

5 

20 

15 

Establish 

702 

^1 

5 

20 

30 

Effects  of 

703 


5 


20 


45 


Turn  Angle 
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TABLE  II  (Concluded) 


Test 

Number 

Gross 
W e i ght 
(lbs) 

Velocity 

(Kts) 

Steering 

(Degrees) 

Braked 
{%  Slip) 

Objectives 

726 

^2 

5 

60 

30 

727 

^"2 

5 

60 

45 

728 

^2 

20 

20 

15 

729 

^2 

20 

20 

30 

730 

^2 

20 

20 

45 

731 

20 

40 

15 

732 

^2 

20 

40 

30 

733 

20 

40 

45 

734 

^2 

20 

60 

15 

735 

^2 

20 

60 

30 

736 

^^2 

20 

60 

45 

The  matrix  evolved  was  generated  assuming  that  gross  weight  changes 
would  be  most  difficult  to  achieve  during  testing,  and  braking  and  steering 
the  least  difficult.  As  a variation  of  parameters  the  matrix  is  correct,  but 
operationally  it  would  seem  more  realistic  to  run  all  tests  of  one  gross 
weight  first,  then  modify  it  to  the  other  value.  Hence,  at  one  site,  all 
gross  weight  tests  would  be  conducted  first  for  all  configurations  of 
velocity,  steering  and  braking,  and  then  all  tests  conducted.  The  tests 
would  then  be  reduced  to  two  32  test  series  at  the  one  site.  Thirty-two 
tests  are  possible  without  significant  pretest  preparation  assuming  the 
taxi  area  does  not  need  preparation  between  tests  and  that  the  gross  weight 
does  not  change  significantly  during  the  time  required. 


Measurements  Required 


The  measurements  required  can  be  classified  as  aircraft  data,  landing 
gear  data,  tire  data,  soil  characteristics,  controls  and  response  data. 

The  aircraft  paraimeters  with  computer  related  symbols  are; 


Gross  Weight 

(W) 

Mass  Moment  of  Inertia  about  the  Vertical  Axis 

(I  ) 

zz 

Taxi  Velocity 

<^o^ 

Ground  Level  to  Thrust  Line 

(UE) 

Engine  Spacing 

(EE) 

Ground  Level  to  Forward  c.  g. 

(U) 

Landing  gear  parameters  are; 

Nose  gear  c.  g.  to  forward  c.  g. 

(L) 

Nose  gear  c.g.  to  Average  c.  g. 

(AL) 

Nose  gear  c.  g.  to  aft  c.  g. 

(LL) 

Nose  gear  c.  g.  to  main  c.  g. 

(F) 

Main  gear  c.  g.  to  forward  c.  g. 

(M) 

Main  gear  c.  g.  to  average  c.  g. 

(AVM) 

Spacing  between  main  gear 

(E) 

Twin  spacing  of  nose  gear 

SN 

Twin  spacing  of  main  gear 

SM 

Tandem  spacing  of  main  tires 

b 'M 

Tire  Data 

Nose  tire  type 

Nose  tire  deflection 

DE 

Nose  tire  diameter 

DN 

Nose  tire  flange  diameter 

DFN 

Total  nximber  of  nose  tires 

NN 

^asmissmsi^sssm 


ag^asssaaKjM 


Main  tire  type 

Main  tire  deflection 

Main  tire  diameter 

Main  tire  flange  diameter 

Total  number  of  n.ain  tires 

Number  of  tire  per  side 

Number  of  tires  in  tandem 

Number  of  tires  in  twin  configuration 


DEM 


DFM 


Soil  Data 


Soil  classification 

Soil  strength  by  cone  index  reading 


NTYP 


Controls  Data 


Braked  slip  as  a function  of  time 

Nose  wheel  steering  angle  as  a function  of  time 


BETA 


Response  Data 

Nose  wheel  drag 
Nose  wheel  side  force 
Nose  wheel  torque 
Nose  wheel  steering  angle 
Nose  wheel  vertical  force 
Nose  wheel  rut  depth 
Left  wheel  drag 
Left  wheel  side  force 
Left  wheel  torque 
Left  wheel  vertical  force 
Left  wheel  rut  depth 
Right  wheel  drag 


TMN 


RML 

SFML 

TMML 

PMR 


RMR 


i 


Right  wheel  side  force 

SFMR 

Right  wheel  torque 

TMMR 

Right  wheel  vertical  force 

PML 

Right  wheel  rut  depth 

ZR2 

Vehicle  side  acceleration  versus  time 

ACVS 

Vehicle  thrust  versus  time 
Path  over  the  ground  versus  time 

THRT 

R 


I 

I 


\ 
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All  of  the  above  are  in  pounds  and  inches  except  for  velocity  which  is 
currently  in  knots  and  side  acceleration  which  is  in  "g”  units. 

Some  of  the  desired  measurements  will  come  from  manufacturers 
data  and  can  be  compiled  before  any  tests  occur.  The  tests  which  will  have 
to  be  measured  by  transducer  or  direct  measurement  are: 

Gross  weight 
Taxi  velocity 
Soils  data 
Controls  data 
All  response  data 

Specific  Measurement  Requirements 

The  following  items  are  discussed  in  order  to  establish  the  range 
of  values  that  can  be  anticipated  and  possible  means  of  collecting  the 
data. 


Gross  Weight 

Gross  weight  has  to  be  found  whether  as  a particular  number  indicative 
of  calculated  values  from  a weights  and  balance  handbook,  or  as  a summation 
of  vertical  forces  measured  at  each  landing  gear  mount.  The  gross  weight 
can  be  obtained  from  strain  gage  data  on  the  struts,  pressures  within  the 
struts,  or  however  it  can  be  easily  measured. 
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Soils  Data 

It  is  desirable  to  have  tests  run  on  clay  type  soil  and  granular  type 

soil..  The  strength  is  determined  by  cone  penetration  tests  and  it  could  be 

that  the  strength  may  dictate  the  gross  weight  ranges  to  be  used.  The 

nomographof  Figure  39  will  be  used  to  establish  the  combination  of  gross 

weight  and  soil  strength  which  is  reasonable.  Standard  bulk  soil  sample 

tests  for  each  .site  are  necessary  for  classification  purposes.  For  the 

range  of  steering  angles  and  velocities  assumed,  it  is  realistic  to  require 

cone  index  data  to  be  collected  at  approximately  250  foot  intervals.  This 

o 

would  generate  about  4 readings  along  the  path  for  the  20  knot,  20  turn. 

Taxi  Velocity 

The  velocity  of  the  vehicle  over  the  ground  should  be  established  for 
the  entire  test  as  a function  of  time.  There  are  several  means  of  achieving 
this  available.  One  would  be  direct  measurement  of  wheel  motion  but  this 
would  be  subject  to  the  problems  of  inherent  slip,  wheel  differential  in 
turns,  surface  roughness  and  dynamic  response.  Another  means  would  be 
optical  tracking  as  by  a photo  theodolite.  This  means  suffers  from  having  to 
differentiate  displacement  measurements  as  well  as  tracking  the  true 
center  of  gravity.  However,  the  optical  tracking  is  suggested  as  the  primary 
means  since  the  path  over  the  ground  is  also  desired. 

Controls  Data 

Braked  wheel  slip  and  nose  wheel  steering  as  functions  of  time  are 
required.  Wheel  slip  specifically,  is  difficult  to  measure.  Consequently 
it  will  be  necessary  to  instead  measure  brake  pressure,  braking  torque 
or  whatever  may  be  available  to  infer  wheel  slip.  Nose  wheel  steering 
should  be  taken  from  a nose  wheel  indicator  if  possible.  It  is  understood 
that  some  steering  devices  have  built-in  limiting  devices  which  do  not  permit 
severe  control  motions.  Therefore,  any  measurements  at  the  control 
wheel  could  be  out-of-phase  and  measurement  should  be  taken  at  the  wheel 
or  mount  as  indicative  of  response  not  control. 
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Response  Data 

All  landing  gear  struts  should  be  instrumented  for  vertical  force, 
drag  force,  side  forces  and  torque.  If  individual  wheel  measurements 
could  be  taken  this  would  also  be  greatly  desirable.  The  specific  force 
levels  for  each  will  have  to  be  based  upon  existing  criteria  limits  unless 
sufficient  data  is  available  ahead  of  time  to  compute  theoretical  values. 

Side  and  drag  force  ratios  of  one-half  can  be  assumed  and  a coefficient 
of  friction  of  0.  8 assumed  to  calculated  generated  torque. 

Rut  depths  should  be  measured  at  each  cone  index  location  and  at 
those  points  where  significant  changes  can  be  observed  due  to  steering 
or  braking  changes.  Because  of  the  many  possible  trajectories  it  will  be 
necessary  to  take  measurements  across  the  ruts  assuming  that  the  terrain 
at  two  to  three  landing  gear  widths  was  unaffected  by  the  passage  of  the 
gear.  Surveying  equipment  will  be  necessary  to  establish  the  location  of 
the  measurements. 

Vehicle  side  acceleration  measurements  can  be  collected  using  an 
accelerometer  with  one  half  *'g"  range  mounted  at  or  near  the  center  of 
gravity  of  the  vehicle  and  aligned  with  the  lateral  axis  of  the  vehicle. 

Longitudinal  measurenr.ents  should  also  be  taken  to  establish  when  an 
equilibrium  condition  is  achieved. 

Engine  thrust  is  requir  data  are  available  to  relate  the  thrust 

to  cockpit  lever  displacemeri-  position  measurements  would  be 

recorded  as  function  of  time.  ' ' is  not  possible,  tachometer  readings 

should  be  recorded  or  engine  * ount  thrust  values  if  possible.,  ; 

j 

i 

} 

( 

I 

f 

i 

I 
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Additional  Comrr.ents 


m-- 


It  is  extrennely  desirable  to  utilize  the  existing  computer  program 
to  study  some  of  the  selected  variations  for  the  particular  vehicle  selected. 

By  doing  that  it  is  possible  to  establish  force  levels  anticipated  and  more 
Importantly  trajectories.  It  is  desired  to  conduct  the  test  for  "first  pass" 
data  collection.  This  infers  that  it  is  possible  to  conduct  tests  that  do  not 
cross  over  the  paths  of  one  another.  Otherwise  surface  preparation  would 
possibly  be  required  after  each  test.  By  examining  the  trajectories  ahead  of 
time  it  may  be  possible  to  establish  not  only  a time  hierarchy  but  a spatial 
hierarchy. 

Test  Facility 

Where  the  tests  can  be  conducted  is  not  known  at  the  present  time. 

Data  Reduction 

It  is  assumed  that  the  majority  of  the  response  and  controls  data 
would  be  collected  using  on-board  recorders  generating  hard  copy  oscillogreim 
records.  The  number  of  channels  required  indicate  probably  two  recorder 
time  synchronized  to  provide  a valid  time  base. 

The  data  channels  can  be  manually  examined  to  determine  where 
equilibrium  conditions  exists  as  seen  in  longitudinal  acceleration,  lateral 
acceleration  drag  force,  side  force  and  vehicle  velocity.  These  locations 
are  then  reduced  using  the  proper  calibration  factors  and  tabulated  in  the 
same  manner  as  shown  on  the  computer  printout  shown  in  Figure  57.  When 
this  type  of  display  it  is  then  possible  to  return  to  the  equations  contained 
within  the  computer  program  and  determine  how  computed  forces  and  sinkages 
compare  with  measured.  Ideally  it  would  be  desirable  to  have  the  data 
reduction  conducted  by  a data  acquisition  system  similar  to  that  used  at 
AEWES  for  tire  tests.  The  system  uses  a computer  to  establish  the  mean 
values  of  the  parameters  and  print  them  out  for  "equilibrium"  conditions. 


AIRCRAFT  PARAMETERS 


Gross  Weight 
Initial  Ground  Vel. 


GW  = 133000.  0 Lbs. 
VO  = 5.0  Knots 


TIRE  DATA  -- 


NOSE 


Tire  Deflection  (%) 

Tire  Diameter  (In. ) 
Flange  Diam.  (In. ) 
Total  No.  of  Tires 
No.  of  Tires/Bogie 
No,  of  Tires--Tandem 
No.  of  Tires  --Twin 
Tire  Type 


DE  = 
DN  = 
DFN  = 
NN  = 


35.00 

38.00 
18.50 

2 


Type  3 


TIRE  AND  C.G.  DISTANCES  --  (IN  INCHES) 


to  FWD  C.G. 
to  AVE  C.  G. 
to  AFT  C.G. 
to  Main  Gear  C.G. 
to  FWD  C.  G. 
to  AVE  C.  G. 


Nose  Gear  C.  G. 

Nose  Gear  C.  G. 

Nose  Gear  C.  G. 

Nose  Gear  C.  G. 

Main  Gear  C.  G. 

Main  Gear  C.  G. 

Ground  Level  to  FWD  C.G. 

Ground  Level  to  Engine 
Between  C.G.  's  of  Left  & Right  Main 
Twin  Spacing  of  Nose  Tires 
Twin  Spacing  of  Main  Tires 
Tandem  Spacing  of  Main  Tires 


SOIL  PARAMETERS 


Cone  Index 
Soil  Type 


Cl  = 192.00 
NTYP=  1 


ENGINE  PARAMETER 

Thrust  Differential  PE  = 


1.00 


MAIN 


DEM  = 
DM  = 
DFM  = 
NM  = 
N1  = 
NMl  = 
NNl  = 
Type  3 


35.00 
55.20 

24.00 
4 
2 
2 
0 


L = 

346.  30 

AL  = 

358.80 

LL  = 

371.30 

F = 

388.00 

M = 

41.70 

AVM  = 

29.20 

U = 

150.00 

UE  = 

153.00 

E = 

172.00 

SN  = 

22.00 

SM  = 

0.00 

SNM  = 

60.00 

Figure  57. 


Master  Program  Input  Requirements. 
(Typed  Facsimile  for  Clarity. ) 
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Should  this  not  be  possible,  manual  techniques  can  be  used  since  there  would 
only  be  64  tests  at  one  site. 

The  reduction  of  the  phototheodolite  data  cannot  be  discussed  until 
more  is  known  about  the  type  of  system  available  (if  that  system  is  used). 

Technical  Order  and  Specification  Developnrient 

As  mentioned  previously,  it  is  very  advantageous  to  conduct  prelinui- 
nary  analyses  for  the  assumed  vehicle  configuration  for  several  of  the  more 
severe  maneuvers.  These  may  indicate  that  force  ratio  will  be  excessive 
and  that  the  braking,  steering  or  velocity  should  be  reduced.  This  in  effect 
is  generating  a set  of  criteria  for  a technical  order  which  will  be  verified 
or  modified  by  test.  Similarly,  *he  forces  developed  may  indicate  that  the 
specifications  currently  used  can  be  modified  in  that  the  force  ratios  for 
wheel  assemblies  or  the  total  vehicle  are  not  realistic  unless  additional 
parameters  are  included  in  the  specification..  That  is,  the  current  specification 
does  not  mention  turning  velocity,  yet  velocity  is  an  important  parameter 
if  side  force  is  to  be  limited  while  conducting  a minimum  radius  turn. 
Additionally,  even  if  the  turn  can  be  accomplished,  the  ruts  developed  may 
not  be  acceptable  because  of  the  required  mission  of  the  unprepared  airfield. 

Preliminary  technical  orders  and  military  specification  modifications 
can  be  evolved  prior  to  testing  and  updated  after  the  test  results  have  been 
evaluated. 
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SECTION  VIII 

MASTER  PREDICTIVE  COMPUTER  PROGRAM 


The  Master  Predictive  Computer  Program  is  composed  of  four 
previously  developed  computer  programs;  TAKEOFF,  LANDING  (TAKE- 
OFF modified  for  landing  aircraft),  TURNING,  and  DPEAC  (aircraift  passes 
capability).  Any  one  of  the  programs  or  emy  combination  of  the  four  pro- 
grams can  be  run  for  a given  set  of  data,  euid  as  many  data  sets  as  are 
required  can  be  used.  A data  set  consists  of  the  values  of  the  program  option 
variables  and  the  input  values  necessary  for  the  particular  program(s)  to 
run. 

The  program  option  variables  are  lOPTl,  IOPT2,  IOPT3,  and  IOPT4 
which  correspond  respectively  to  the  TAKEOFF,  LANDING,  TURNING, 
and  DPEAC  programs.  By  setting  the  appropriate  option  variable  equal  to 
an  integer  1 (one)  for  TAKEOFF,  TURNING,  and  DPEAC  and  equal  to  an 
integer  2 (two)  for  LANDING  the  designated  program(s)  will  be  executed. 

If  a program  is  not  to  be  run,  an  integer  other  than  those  mentioned  above 
must  be  supplied  as  input, 

DASET,  the  number  of  data  sets  to  be  run,  is  the  first  variable  read. 
Next  to  the  program  option  variables  (lOPTl,  IOPT2,  IOPT3,  IOPT4) 
are  read  followed  by  the  variables  common  to  ail  four  programs.  The  rest 
of  the  variables  read  in  are  specific  to  particular  programs  as  illustrated 
in  Figure  58. 

All  of  the  variables  used  in  the  master  program  have  been  previously 
used  and  defined  in  the  TAKEOFF,  TURNING,  and  DPEAC  programs 
except  for  the  following: 

THMAX  - moiximum  aircraft  thrust 
VOL  - initial  velocity  for  LANDING 

VTOL  - dummy  takeoff  velocity  variable  in  LANDING  (VTOL  must  have  a 
greater  value  than  VOL  for  LANDING  to  run) 
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AOL  - initial  acceleration  for  LANDING 

SOL  - initial  rvinway  distance  for  IjANDING 

V03  - initial  velocity  for  TURNING  (VO  in  original  TURNING) 

DT3  - time  increment  for  TURNING  (DT  in  original  TURNING) 

An  example  of  the  data  cards  required  to  run  the  master  predictive 
program  is  shown  in  Figure  59  for  the  C-130E  aircraft.  The  program 
option  variables  in  this  example  (data  card  number  Z)  have  been  supplied 
with  values  to  run  all  four  programs.  As  an  illustration  of  using  combinations 
of  the  programs,  the  program  option  variables  data  card  shown  beneath 
the  data  set  has  values  which  will  run  only  TAKEOFF  and  LANDING. 

The  outputs  generated  by  the  selected  routines  are  identical  to  those 
previously  discussed.  That  is,  the  turning  portion  provides  output  as  seen 
in  Figure  56,  and  the  number  of  passes  are  presented  in  Reference  ZO. 

The  main  purpose  of  developing  the  master  program  was  to  generate  one 
program  which  would  use  the  saime  variables  in  essentially  three  different 
digital  routines.  This  aspect  was  accomplished  but  the  outputs  were  not 
modified  to  have  a "standardized”  appearance. 
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Figure  59.  Input  Data  Cards  for  Master  Program. 

(Typed  Facsimile  for  Claritv  ) 
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SECTION  IX 


APPLICATION  TO  CURRENT  SPECIFICATIONS 


A.  TURNING 

Military  Specification  MIL-A-008862A  , Airplane  Strength  and 
Rigidity,  Landing  and  Ground  Handling  Loads,  contains  the  definitioii  of 
a turn  and  formula  to  be  used  to  establish  turning  loads.  Under  paragraph 
3.3.2  Turning  the  turn  is  to  be  achieved  by; 

a)  Unsymmelrical  thrust  or  nose  wheel  steering 

b)  Unsymmetrical  thrust  or  nose  wheel  steering  with  symmetrical 
braking,  and 

c)  Differential  braking. 

The  parauneters  that  generate  a turn  are  the  braking,  differential  thrust  or 
nose  wheel  steering.  If  we  tabulate  the  various  values  these  parameters 
can  have,  there  are  twelve  combinations.  For  example 


Brakes 

TTirust  { 

Symmetrical) 

Thrust  (Unsymmetrical) 

Trail 

Steering 

Trail 

Steering 

None 

1 

3 

2 

4 

Symmetrical 

5 

7 

6 

8 

Differential 

9 

11 

10 

12 

Cells  1 and  5 are  trivial  combinations  since  either  no  brakes  or  symmetrical 
brakes  with  a trailing  nose  wheel  and  symmetrical  thrust,will  not  turn 
the  vehicle. 

If  we  consider  turns  achieved  by  only  one  control  parameter,  there 
are  three  possibilities.  Cells  2,  3,  and  9 are  for  a turn  generated  by 
thrust  alone,  steering  alone,  and  differential  braking  respectively. 


Two  control  parameter  turns  are  those  of  cells,  4,  6,  7,  10,  and  11. 

These  are  asymmetrical  thrust  with  nose  wheel  steering  and  symmetrical 
brakes,  nose  wheel  steering  with  symmetrical  brakes,  differential  braking 
with  asymmetrical  thrust  and  with  nose  wheel  steering.  Symmetrical 
brakes  introduces  a unique  turning  condition  because  they  influence  the 
vertical  forces  at  the  landing  gear  by  changing  the  pitching  moment. 
Therefore,  even  though  symmetrical  braking  does  not  change  the  direction  of 
the  aircraft,  its  combination  with  other  parameters  changes  the  turning 
response.  All  three  parameters  are  used  in  cells  8 and  12  where  differential 
thrust,  nose  wheel  steering,  and  symmetrical  and  differential  braking  are 
possible. 

Current  specifications,  as  listed  above,  consider  only  cells  2,  3,  6, 

7,  and  9.  Combining  these  with  the  trivial  cells  and  eliminating  these  from 
the  total,  there  are  five  conditions  not  covered  by  the  current  specifications. 
These  are: 

Asymmetrical  thrust  with  steering 

Asymmetrical  thrust  with  steering  and  symmetrical  brakes 

Asymmetrical  thrust  with  unsymmetrical  brakes 

Differential  braking  with  steering 

Differential  braking  with  steering  and  unsymmetrical  thrust 

From  the  table  it  is  apparent  that  those  not  currently  specified  are  those 
of  the  most  extreme  maneuvering  inputs,  those  with  the  greater  ntimber 
of  control  parcimeters.  From  very  limited  data,  it  appears  that  these 
conditions  would  be  those  that  would  create  the  smallest  operational 
turning  radius  and  would  therefore  be  very  desirable  maneuvers  for 
aircraft  operating  in  forward  areas.  The  data  presented  in  Section  III 
indicate  that  for  a particular  set  of  input  conditions  to  a C-130E,  the 
minimum  field  width  could  be  reduced  by  20  percent  by  using  differential 
braking  with  nosewheel  steering.  The  thrust  required  to  overcome  the 


braked  drag  is  nearly  all  of  the  thrust  available  and  therefore  in  this 
case  differential  thrust  would  not  even  be  possible.  Therefore,  it  appears 
that  while  some  of  tlie  combinations  are  at  least  theoretically  possible, 
they  may  not  in  fact  be  practically  possible.  Nonetheless,  there  are 
additional  turning  maneuver  operations  which  could  conceivably  be 
considered  for  inclusion  into  the  specification. 

Under  Paragraph  4.  6 Turning,  of  the  standard,  several  formulas  are 
given  for  calculating  the  loads  for  outside  and  inside  gear  based  upon  ver- 
tical load  factor,  N , side  load  factor,  N , and  weight  distribution  based 

6 

upon  the  geometry  of  the  landing  gear.  This  is  consistent  with  assuming 
that  the  vehicle  is  in  a steady  state  turn  where  side  forces  balance  the 
inertial  response  and  the  vertical  forces  are  generated  by  inertial  response 
cind  rolling  moment. 

From  previous  discussions,  it  is  apparent  that  rolling  drag  can  become 
very  large.  Therefore,  it  is  not  necessarily  realistic  to  ignore  the  effects 
of  rolling  drag  developed  in  soils.  The  specification  indicates  that  for  bare 
soil  fields,  each  unbraked  wheel  shall  have  a drag  reaction  of  0.  2 times 
the  vertical  reaction.  If  this  is  observed,  then  it  is  only  reasonable  to 
have  a forward  load  factor,  N , which  would  be  used  to  calculate  an 

Li 

additional  contribution  to  the  vertical  forces  on  each  landing  gear.  The 
current  specified  value  of  0.8  times  the  vertical  reaction  for  lateral  load 
(drag  and  side  force  resultant)  seems  reasonable.  The  drag  ratio  at  that 
level  would  add  significantly  to  the  nose  wheel  loads. 

The  specified  side  load  ratio  of  0.  5 for  both  vehicle  and  any  wheel 
seems  reasonable  for  vehicle  and  main  gear  based  upon  very  limited  data. 
The  C130E  example  indicates  that  the  vehicle  side  force  ratio  C3in  only  be 
approached  at  high  taxi  speeds  where  the  turning  radius  becomes  very  large. 
Hence,  considering  the  desire  to  keep  taxi  speeds  down  and  turn  in  a 
reasonable  space,  the  vehicle  side  force  ratio  is  probably  consistent  with 
other  turn  parameters. 


The  side  force  ratio  as  applied  to  any  one  wheel  may  be  more  suspect. 
The  first  two  specified  means  of  turning  are  with  nose  wheel  steering.  Test 
data  indicates  that  the  limiting  value  may  be  exceeded  depending  upon  the 
turn  angle  and  sinkage.  Additionally,  the  peak  side  force  on  the  nose  wheel 
may  occur  because  of  the  steering  rate,  not  the  steering  angle.  Therefore, 
it  appears  that  the  nose  wheel  may  require  a special  analysis  to  predict 
its  maximum  loading  condition.  Either  that  or  retaining  the  present  maxi- 
mum side  force  ratio  will  require  that  steering  rate  and  forward  velocity  be 
examined  to  calculate  their  effects.  This  is  necessary  since  the  turn  angle 
is  dictated  both  by  the  steering  rate  and  the  forward  velocity. 

Present  turning  specifications  do  not  reflect  all  possible  configurations 
of  turning  mechanisms  possible.  The  equations  available  do  provide  real- 
istic load  factors  for  vehicle  side  force,  but  do  not  reflect  the  pitching 
moment  due  to  drag  and  thrust.  Nosewheel  steering  can  provide  unique 
loading  conditions  due  to  steering  rate  and  forward  velocity  auid  these  factors 
should  be  mentioned.  Additionally,  with  current  equations  and  nomograms 
available  for  the  solution  of  sinkage  and  drag,  it  is  possible  to  obtain  better 
estimates  of  rolling  and  braked  drag  rather  than  assuming  two  distinct 
values.  Lastly,  because  of  the  effects  of  speed,  not  only  for  drag  but  upon 
turning  rate,  velocity  should  be  introduced  in  order  to  call  attention  to  its 
presence. 

B.  MULTIPASS 

The  Design  Handbook,  AFSC  DM-21,  under  Design  Notes  4C2,  Un- 
surfaced Airfields,  lists  a method  to  be  used  in  evaluating  relative  ground 
flotation  performance.  The  procedure  enables  calculation  of  the  permissable 
number  of  passes  an  aircraft  can  make  over  heavy  clay  soil  prior  to  the 
formation  of  a three -inch  rut.  The  procedure  is  restricted  to  one  soil  type 
and  one  soil  failure  criteria.  Data  collected  from  this  program  indicate 
that  rut  depth  cein  be  predicted  as  a function  of  soil  type  and  strength,  and 
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number  of  passes.  The  data  are  currently  restricted  to  a limited  number  of 
pas  ses,  but  it  is  indicative  that  the  depth  can  be  predicted  since  the  three- 
inch  rut  may  not  be  limiting  for  large  diameter  tires. 

Siiioe  the  prediction  of  rut  depth  relies  upon  first  pass  sinkage  and  rut, 
it  is  possible  that  additional  nomograms  or  equations  could  be  included  in 
the  multipass  analysis  currently  used.  The  technique  requires  the  calculation 
of  tire  contact  area  and  single  wheel  load.  The  solution  of  first  pass  sinkage 
and  rut  depth  then  only  requires  knowledge  of  soil  type  and  strength.  With 
this  as  a starting  point,  rut  depth  for  a given  number  of  passes  in  the  same 
path,  or  adjacent  paths  could  be  determined.  In  this  manner,  it  is  possible 
to  determine  an  allowable  number  of  passes  for  a variable  rut  depth.  If  it 
can  be  shown  that  the  permissable  rut  depth  is  dictated  by  limiting  drag  or 
structural  interference,  then  the  number  of  passes  would  be  dictated  by 
tire  size  and  deflection,  soil  properties,  and  the  particular  vehicle  geometry, 
not  a selected  rut  depth. 

A revised  means  of  calculating  allowable  passec  would  be: 

(1)  Establish  msotimum  permissable  rut  depth  based  upon  landing 
gear  physical  configuration,  external  stores,  or  vehicle 
consideration  such  as  wheel  well  doors. 

(2)  Establish  mziximum  permissable  rut  depth  based  upon 
structural  strength  restrictions  or  drag  restriction  determined 
by  thrust  available. 

(3)  Select  maximum  permissable  rut  depth  from  the  above  information. 

(4)  For  the  specific  landing  gear,  tire,  soil  and  vehicle  weight 
distribution,  calculate  the  number  of  passes  possible.  This 
can  assume  multiple  passes  in  the  same  lanes  or  a specific 
distribution  and  time  phasing  between  passes.  Braking  could 
be  included. 


In  siunmary,  multipass  design  criteria  can  be  updated  so  that  the 
appropriate  specifications  would  require  the  allowable  number  of  passes 
to  be  found  based  upon  new  predictive  equations.  These  would  reflect  soil 
strength  and  type,  and  the  time  history  of  the  passes  made. 

C.  TAKEOFF  AND  LANDING 

Takeoff  and  landing  performance  on  unpiepared  runways  is  specified 
in  terms  of  rolling  or  braked  drag.  Rolling  drag  is  assumed  to  be  0.  2 times 
the  vertical  reaction  while  braked  drag  is  0.8  times  the  vertical  reaction. 
These  can  be  altered  to  easily  reflect  all  influencing  parameters  except 
forward  velocity.  There  are  nomograms  which  permit  calculation  of  rolling 
drag  as  a function  of  tire,  soil  and  single  wheel  load.  Hence,  the  rolling 
drag  coefficient  for  the  particular  vehicle  could  be  evaluated  using  specified 
procedures,  rather  than  having  specified  values  for  all  aircraft.  Similarly, 
the  results  of  the  rolling  drag  analysis  can  be  used  in  braked  drag  equations 
to  determine  the  effect  of  braking  upon  landing  rollout. 

Of  particular  significance  in  calculating  braked  rollout  would  be  the 
evaluation  of  braked  rut  depth.  The  braked  rut  is  greater  than  the  free 
rolling  and  the  performance  of  multiple  landings  and  takeoffs  could  be  severely 
restricted  unless  this  is  recognized.  Consequently,  any  landing  and  take- 
off specification  should  reflect  not  only  the  drag  ratio  analysis,  but  rut 
depth  generation  as  it  influences  multipass  operations. 
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SECTION  X 

CONCLUSIONS  AND  RECOMMENDATIONS 
A.  CONCLUSIONS 

The  reported  results  of  the  past  years  research  have  generated  the 
following  conclusions. 

Turned  Tire  Phenomena 

1.  Drag,  sinkage  and  rut  depth  increase  with  turn  angle.  The 
increases  are  greater  in  cohesionless  soils  than  in  cohesive 
soils  for  a given  load  and  soil  strength. 

2.  The  lateral  load  ratio  can  amd  in  some  cases  will  exceed  one -half 
where  the  turn  angle  becomes  large, 

3.  Pneumatic  trail  is  negative  (behind  the  tire  axis)  and  decreases 
with  turn  angle  in  cohesive  soil.  Pneumatic  trail  is  positive  and 
increases  with  turn  angle  for  cohesionless  soil.  Therefore, 

the  trail  generates  a stable  condition  in  cohesive  soil  and  un- 
stable in  the  cohesionless. 

4.  Realistic  results  are  generated  by  an  analytical  model  of  a 
turning  aircraft  using  predictive  equations  from  tire  test  data. 

5.  Computed  turning  results  indicate  that  current  specified  nose 
wheel  steering  angles  generate  turn  angles  greater  than  those 
used  in  the  laboratoiry  tests.  Additionally,  load  strength  ratio 
limits  were  exceeded. 

6.  A limited  parametric  study  indicates  that  the  optimum  means  of 
turning,  considering  force  ratios  only,  is  nose  wheel  steering 
only,  nose  wheel  steering  with  engine  thrust  differential,  and 
differential  braking  last.  If  am  extremely  tight  turn  is  required, 
only  an  analysis  considering  the  vehicle  configuration  can 
determine  the  best  balance  of  steering,  braking,  and  thrust. 


Multtpaf  Rcaponse 


1.  Rut  depth  generated  by  a single  rolling  tire  passing  over  the 
same  path  in  cohesive  soil,  can  be  calculated  by  a linear  change 
from  the  first  pass  rut  depth  for  approximately  twenty  passes  not 
exceeding  a three -inch  cumulative  rut  depth, 

2.  Drag  forces  measured  on  the  same  single  free  rolling  wheel 
passing  over  the  same  path  do  not  significantly  increase  in  a 
cohesive  soil. 

3.  A single  wheel  which  alternates  path  and  braking  generates  a rut 
depth  which  can  be  calculated  using  superposition  of  constant 
coefficient  effects  all  related  to  first  pass  sinkage  and  rut  depth. 

Start-Up  Forces 

1,  The  limited  data  examined  indicate  that  the  start-up  drag  force 
ratio  can  be  approximated  as  a function  of  the  first  pass  free 
rolling  drag  ratio. 

2.  The  data  indicate  that  there  is  a trend  which  follows  the  load 
strength  ratio,  a /Cl,  but  there  are  other  effects  more  directly 
related  to  total  force  and  cone  index  alone  which  are  not  yet 
defined. 

Roughness  Effects 

It  is  possible  to  model  the  soft  tire /soft  soil  response  using  lumped 
parameter  elements  as  well  as  soil  properties.  The  computational  time 
may  be  excessive  for  a "real  world"  simulation. 

B.  RECOMMENDATIONS 
Turned  Tire  Phenomena 


1.  Additional  tests  are  desired  for  tires  operating  at  greater  turn 
angles  with  and  without  braking.  It  is  not  known  how  the  braked 


drag  is  influenced  by  the  siniage  created  by  the  turn  amgle.  This 
effect  can  be  significant  as  shown  by  the  computed  results. 

2.  Full-scale  testing  of  an  instrumented  vehicle  as  developed  in  the 
Test  Plan  of  Section  VII  is  necessary  to  validate  the  digital 
program  developed  from  tire  test  results. 

Multipass  Response 

1.  All  results  generated  during  this  research  e'fort  were  for 
cohesive  soils.  Therefore,  the  next  step  required  is  that  of 
conducting  another  test  series  using  a cohesionless  soil  for 
similar  test  parameter  ranges. 

2.  The  direct  and  indirect  rolling  and  braking,  rut  depth  coefficients 
exhibit  an  apparent  reversed  at  a particular  value  of  P/CI,  Tests 
should  be  conducted  over  that  remge  to  validate  the  limited  data 
available. 

Start-Up  Forces 

1.  A test  series  should  be  conducted  to  isolate  the  effects  of  time, 
tire  pressure,  load  strength  ratio,  emd  total  force  upon  start-up 
force. 

Roughness  Effects 

1.  The  lumped  parameter  values  for  the  computer  model  should  be 
established  for  several  tire  configurations  using  existing  test 
data.  There  are  test  results  available  to  establish  tire  coefficients, 
rolling  tires  on  rigid  surfaces,  and  to  establish  soil  coefficients, 
rigid  tires  on  soft  soils.  After  these  have  been  developed,  selected 
"roughness"  configurations  can  be  studied. 
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AIRCRAFT  PARAMETERS 


GROSS  WEIGHT 

GW  = 

133000. 

0 LBS. 

INITIAL  GROUND  VEL. 

VO  = 

c 

.0  KNOTS 

TIRE  DATA  — 

NOSE 

MAIN 

TIRE  DEFLECTION n:) 

DE  = 

35. 

00 

OEM  = 

35.00 

TIRE  DIAMETER  (IN.) 

ON  = 

38. 

00 

DM  = 

55.20 

FLANGE  OIAM.  (IN.) 

OFN  = 

18. 

50 

OFM  = 

24.00 

TOTAL  NO.  OF  TIRES 

NN  = 

2 

NM  = 

4 

NO.  OF  TIRES/BOGIE 

N1  s 

2 

NO.  OF  TIRES— TANDEM 

NMl  = 

2 

NO,  OF  TIRES— TWIN 

NNl  s 

0 

TIRE  TYPE 

TYPE  3 

TYPE  3 

TIRE  AMO  C.G.  DISTANCES— (IN  INCHES) 


NOSE  GEAR  C.G. 
NOSE  GEAR  C.G. 
NOSE  GEAR  C.G. 
NOSE  GEmR  C.G. 
MAIN  GEAR  C.G. 
MAIN  GEAR  C.G. 
GROUND  LEVEL 
GROUND  LEVEL 
BETWEEN  C.G.*S 
TWIN  SPACING 
TWIN  SPACING 
TANDEM  SPACING 
SOIL  PARAMETERS 
CONE  INDEX 
SOIL  TYPE 


TO  FWO  C.G. 

L 

s 

346, 30 

TO 

AVE  C.G. 

AL 

s 

358.80 

TO 

AFT  C.G. 

LL 

r 

371.30 

TO 

MAIN  GEAR  C.G. 

F 

s 

388.00 

TO 

FWO  C.G. 

M 

s 

41.70 

TO 

AVE  C.G. 

AVM 

r 

29.20 

TO 

FWO  C.G. 

U 

= 

150.00 

TO 

ENGINE 

UE 

5 

153.00 

Of 

LEFT  L RIGHT  MAIN 

E 

= 

172.00 

OF 

NOSE  TIRES 

SN 

= 

22.00 

OF 

MAIN  TIRES 

SM 

0.00 

OF 

MAIN  TIRES 

SNM 

= 

60.00 

CI=  192.00 
NTYP=  1 


ENGINE  PARAMETER 

THRUST  DIFFERENTIAL  PE= 


1.00 


C-130  Configuration  Vehicle  Characteristics 


fllRCPftPT  PAPAMETERS 


GROSS  WEIGHT 

GW  = 310!: 

LBS 

' • 

INITIAL  GPOUNO 

VEL 

. VO  = 5. 

C KNOTS 

TIRE  OATA  — 

NOSE 

MAIN 

TIRE  0£PLECTION(X» 

OE  = 32.0- 

OEM 

- 

32. CO 

TIRE  niAMETER 

(IN. 

) ON  = 23.73 

DM 

S 

34.55 

flange  oiam. 

(IN. 

) OFN  = 11.81 

OFM 

= 

18.75 

TOTAL  NO,  OF  TIRES 

NN  = 1 

NM 

s 

2 

NO.  OP  TIRES/BOGIE 

N1 

s 

1 

NO.  OF  TIRES  — 

TANOEM 

NMl 

C 

NO,  OF  TIRES  — 

TWIN 

NNl 

r 

G 

TIPF  TYPE 

TYPE  7 

TYPE 

7 

TIRE  AND  C.G.  OISTANCFS— (IN  INCHES) 

NOSE  GEAp  C.G. 

TO 

FWD  C.G. 

L 

= 

184.70 

NOSE  GEAR  C.G. 

TO 

aW  C.G. 

AL 

s 

188,25 

NOSE  GEAR  C.G. 

TO 

aft  C.G. 

LL 

= 

191.79 

NOSE  GEAR  C.G. 

TO 

MAIN  GEAR  C.G. 

F 

212,79  ■ 

MAIN  GEAR  C.G. 

TO 

FWD  C.G. 

M 

= 

21.i:i 

main  gear  C.G. 

TO 

AVE  C.G, 

AVM 

24.55 

GROUNO  LEVEL 

TO 

FWO  C.G. 

U 

= 

89,10 

GPOUNP  LEVEL 

TO 

ENGINE 

uE 

= 

113T43 

8ETWEFN  C.G. *S 

OF 

LEFT  1 RIGHT  MAIM 

E 

206.85 

TWIN  SPACING 

OF 

NOSE  tires  ' 

SN 

= 

o • 0 Q 

TWIN  SPACING 

OF 

MAIN  TIRES 

SM 

= 

-.CO 

tanoem  spacing 

OF 

MAIN  TIRES 

SNM 

j.OQ 

soil  parameters 

CONE  INDEX 

CT=  19'2.Cf? 

SOIL  TYPE 

NTYP=  1 

- - - - - 

fmqine  parameter 

- - 

- -■ 

- - - - 

THRUST  DIFFER" 

NTIAL  Pc=  1.:!! 

Attack  Vehicle  Characteristics 
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AIRCRAFT  PARAH£TFRS 

GROSS.  MElGHr  . GH  .=  . _2±aCX  U G.  J_3iS^ ..  ..  _ . 

INITIAL  GROUNO  V£L.  VD  = 5.3  KNOTS 

TIRE  DATA—  NOSZ  MAIN 


TIRE  JlEFLEfll ION (X.)  ...  _DE_=. 32. 3i_..HEM„=...  .32.3  0 


TIRE  DIAMETER  (IN.)  ON  = 
FLANGc  LIAM.  (IN.)  OFN  = 
TOTAL  NO.  OF  TIRES  NN  - 
NO.  OF  TIRES/aOaiE 
NO.  OF  TIRES  — TANOEM 

.JiO..  jaf-J.XRESr-.miN  . - 

TTRE  TYPE  TYPF^^ 


‘‘O.cc 

19.50 

2 


DM  = 
OFM  = 
NM  = 
N1  = 
NMl  = 
...NNl  = 
TYPE 


40.00 

19.50 

6 

4 

4 

. . 4- 


TIRE  ANU  C.G.  DISTANCES— (IN  INDIES) 


..NOSE.  GEAR-  C RHD  C.G^_. 

. L- 

S 

-.42J-.aO 

NOSE  GEAR  C.G,  TO 

AVE  C.G. 

AL 

= 

44, ,20 

NOSE  GEAR  C.G.  ..TO 

AFT  C.G.  . .. 

- LL 

= 

-452.60 

NOSE  GEAR  C.G.  TO 

MAIN  GEAR  C.G. 

F 

S 

487.00 

.MAIN  -GEAR  C_G.  TO 

FHD  C.G* 

M 

59.20 

MAIN  GEAR  C-G,  TO 

AVE  C.G. 

AVM 

= 

46.80 

GRQUNQ  LEVEL  TO 

-5)(D  C..£l.  _ - - . 

-U . 

17o,a.a 

GROUND  LEVEL  TO 

ENGINE 

UE 

241.00 

BETHEEN  C.G,  »S.OF 

LEFT  & RIGHT  MAIN 

r 

= 

223.00 

THIN  SPACING  OF 

NOSE  TIREi- 

SN 

s 

31.00 

THIN  SPACING  OF 

MAIN  TIRES-  - - 

s:i 

60.00 

TANDEM  SPACING  OF 

MAIN  TIRES 

SNM 

35.00 

CONE  INDEX 

CI=  192.00 

SOIL.  TYPE 

NTYP=  1 

_ 

ENGINE  PARAMETER 
tNEUSJ_OI£££  8EN1IAL  _E£=_.. 


L.JiX- 


Advanced  Tvanspcrt  Characteristics 


•••Mlh’C  SOUfiCf-LAHCUACE  PROGRAM***. 
VEFSIOM  t 1C/C1/6A  HOO  LGyLL  COOl 


SQFl  TIRE  ROLLING  OH  A SfFI  EURFACt- 
CQN(VO,ZO>P.OP.HBiHP) 

COH(HN,HM,RM.IIO  

CON(«T  .RT.CTtBn 

XON<HS,<S,CS.CSS.HSI 

CON(COHS,IFS» 

calculate  N££0£0  CONSTANTS-  - 


MT 
MM 
HP 

OH 

AZ 


MT/3J6. 

HH/3B6. 

HP/386. 

MS/3J6»- 

VO/PM 


BI*:.1571*RH  

3 CALCULATE  ACCELFRATION  OF  MHEiL  HUB  VERTICALLY 

ZOZMi  (-P*TSUH+FSUHt/<MH*MP)-386.  

31  CALCULATE  SUN  IF  TAHG'NTIAL  FORCES 

TSUM  ISUK1*I3UH2*ISUH3 

TSUHl  Tl*COS(rHll*T2*COS<IH2l*T3*COS«TH3»*Tfc*COS(TH4) 

TSUKZ  Ts’COS  (TH5l*r6*C0S(TH6)  *T7*C05  (TH7»*TB*C0S  (TH6).- 

TSUH3  T9*COStT  '1«T10*COS«TH10) 

32  . CALCULATE  SUN  OF  SPOKE  FORCES  . 

FSUM  FVTl*SIIurMl»*FVT2*SINITH2»*FVT3*SIN  (TH31+FSUM1 

FSUMl FVTL*S  IN  < TH4 » ♦FVT5»SIN  f rM5l*FVI6*SIMttM5*>FSUH2- 

FSUM2  FVT10*SIN(TH10I 

33  TANGENTIAL  FORCE  PcRNITTEO  AFTER  IMPACT 

T1  LSHTGO.U’.O.I 

12  L3N(G02,T2J, J.J_ 

T3  LSM(G03.T3B,a.l 

IL—,. LSHCGOA.TvB.a^l 

T5  LSM(Ga3.TSOi!).l 

T6  LSM(G06,T6B,a,»- 

T7  LSM(G07,T79,a.l 

TB  LSM(C03,T8B.3.» 

T9  LSM(G09,T9B.O.I 

tlJ L3M«GOl,,TlCa,0.l- 
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RELATIVE  VELOCITY  DICTATES  SIGN  OF  T 

- -TIB  FSH(RVEL,'TIA,3.,UAI 

T20  FSH(RV£L2,-T2A,(1.,T2A1 

- -T30  FSH(RVEL3,-T3A,0.,T3A» 

T4B  FSM<RVEL4,-T*tA,S,,T<AI 

X5B FSM(.RVELS,-T5A,a.^TSAl 

T6B  FSHIRV:L6,-T6A,a.«T6AI 

T7B  F3H(RVEL7,-T7A,l!.,T7Al 

TJB  FSMIRV£L8,-TaA,0.,T8A» 

T93  FSHIRVEL9,-T9A,0.,T9AI - 

Tloe 


FSH(RV£L10t-T10A,U.,T10AI 
_ RELAIIVt  VELOCITY  AT -SPOKE  SURFACE- 
RVEL 
.-  RVELE 
RVFL3 
RVEL4 
RVELS 
_JTV£L6  . 

RVtL7 
RVEL8 
RVEL9 
RVELIO 


R*10THI-10XH*SIN(TH1I 

R2*10THl-lDXH»SINnH2> 

R3*10TH1-1')XH*SIN(TH3> 

R4*10TN1-10XH*SIN(TN4* 

R5*10THl-10XH*SimTHE» 
Po*10THl-10XM*SimTH61 
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Rr*10THl-l()XH*SIM(lH7l 

P’>»10TH1-10XH*SINCIH8» 

R9*10THl-inXK*SIMIIH9» 

R1Q*10TH1-10XH*SIN(TM10I  

MAGNITUDE  OF  FRICTIONAL  FORCf  T 
TIA  (AZ*COHS*Nl*IFSI*MUSL  


TJA  «AZ*C0HS*N2»IFSJ»MUSL2 

....  I3A  . _ (AZ*C0M3*NJ*  IfSI ‘HUSLI _ 

TAA  (AZ*C0M5*N-»IFSl*HUSL4 

...  ISA  (AZ*C0HS*N5*IFSI»HUSL5 

T6A  (AZ*COHS»Mij»IFSI ‘MUSL6 

I7A (AZ*CQH3«NZ»IFSI  «HUSLZ 

T»A  (AZ*C0HS*‘I««IFS1»HUSL» 

..  T9A  (AZ*C0H3*N9»IFS1»MUSL9 

TlOA  (AZ’COHS+'Ufl^IFSJ^MUSLlQ 

CU3E  ROOT  OF  SLIP  ..  

HUSL  £lfP(a.313,SLIP> 

“USLi EXP « 0 . 3 33 . SLIP21 

MUSC3  EXPia. 133.SLIP3I 

PUSL4  ..  EX^O.  333.SLIP41 

MUSL5  EXP(3. 533, SLIPS* 

...  . MgSLS  £XP(3. 333, SLIPS! 

MUSL7  EXO<3. 333,SLIP7I 

!<USLa £XP(  3.  333,SL1P«1 

MU3L9  EXP(0. 333.3LIP9I 

•5USLU  iXO(J.  333.SLIPlQ> 

SLIP  IS  RtLATIv;  VELOCITY  OVER  R OHEGA 

SUP  A9S(RVEL/U0TH1»R1 1 

SLIPZ  A9S(RV-:L2/(10IH1»R2H 

SLIP3  ..  AQS(RViL3/(10rHl»R3ll 

SLIP4  AbS(9V£L4/(ialHl»R4ll 

. SLIPS  A3S(9VEL5/510THI»R5)I 

SLIPS  A9S(RV;LS/(iaTHl»R6! I 

. . SLIP7  ABS(9VEL7/(10TH1*R7!» 

SLIPS  AnS(RVEL8/(13THl»R8*» 

SLIP.9  AU3<9V£L9/«1JTH1?JI9AI 

SLiPio  Aos<Rv-;Lia/(iorHi»Rioii 

FORCES  IN  THE  S°OKE  _. ...  

FVTl  <KT*  <7U-ZWU«r,T»(10ZTl-10ZHll)/SIN«IHll 

FVIZ  (KT*  (7T2-ZH2)  *01  *aD7T2-lOZHZ!)/SI“M !,!?)_ 

FVT3  (XT»<7f3-7<3l»CT*(102T3-l0ZM3ll/SIN<IH3) 

F.VI.4 (KT»(7I4-ZM41  »CI*aDZT-,'.10ZM4ll/SIHlIJ34» 

FVT5  (KT*(ZI5-ZM5ltCT*a0ZT5-10ZH5!!/SINIIH5) 

_..  FVT6  (KT*  (7T6-ZM6J*CT»UDZTS-10ZH6II/SIN(THS) 

FVI7  (KT*<ZT7-ZM7l*CT*(10ZT7-10ZH7)>/SI'mH7> 

_..  FVT8  ._  (KT»(ZT8-7W8I  ♦CT»(10ZTe-10ZWSI)/SI!mH9) 

FVT9  (KT»(ZI3-ZH9!*CT»a0ZI9-lOZH9l!/SISITH9) 

FVT15 <KT*  (ZT10-ZH10)»CT».aOZ.TJ.6rlOZH103i/SlN(.TJtLll 

CALCULATE  HUB  RELATIVE  VELOCITY 

SPCKE  VELOCITIES  AFTER  IMPACT 

19ZM1  LSH(GO,lDZHlB,a. ) 

10ZH2  LSH(GO2,19ZHZB,0.) 

10ZK3  LSH(G03,10ZH3a,O.I 

1QZH4 LSM(G04,10ZH49,Oj.l 

lOZNS  LSH(GOS.10ZMSB,O.I 

-.lOZHS  LSH(GOS,10ZHSn,O.I 

10ZH7  LSH(G07,13ZH7B,0.! 

10ZM8  LSW(G09.13ZH89,O.I 

10ZM9  LSH(G09.10ZM9B,3.> 

lOZHlt  L5H(GO10,13ZHAB.Q.). 

INTEGRATE  ACCELERATIONS  AFTER  IMPACT 

lOZMlB  _ l!tT(2CZHl,Hl,TRU£.GOI 

10Z42B  INT(20ZVI1,HZ,TRUE,G02I 

10ZH3B  INTtZDZNliHStTRUE.GOSI . . ..  .. 

I0ZH4P.  IHT  (ZOZH1,V4,TRUE,G04I 

30ZH5B  INT(20ZH1,H5, TRUE, COS) 

lOZWSB  INT(20ZM1,R6, TRUE. COS) 

10ZM7B  INT(20ZV1,M7.TRUE,G07! 

10ZM8B  IMT(20ZM1,H8,TRUE,G08I 

10ZR9B  INT(20ZM1,H9, TRUE, GOB!  

IBZMAB  lNT(Z0ZHl,wn,TRUE,GO10! 

TRACK  ANO  STORE  INERTIAL  VELOCITY 


W1  TASUOZdlA.I'll.a.J 

H2  TASU02KU.IM2.3.J  

>43  TAS(102ulA,t'43,J.I 

Uli  . TASU02M<A.IU‘>,0.i 

m;  tasiiozmia,im$,3.i 

H6 IASI102HU.IN6.Q.3 

U2  TAS(1C2M1A,1N7,3.I 

M3  TASaoZMlA.I^SiO.I  

M9  TAS(102MlA.IN9,a.l 

MIC  TAS(102HlA,I>9ia.O.I 

lOZMlA  IHT(20ZM1.ZO» 

calculate  spoke  OISPLAC£M£NI_AF.I£a.IMPA 

ZMi  TUT  (107>41,9.  ,I9U£,G0I 

ZM2  HT(10ZW2,0.,3RU£.C02» 

ZM3  INru0ZM3.3..TRUE>GO3l 

ZMA  :NT(10Z9A,a.,TRU£.G0l.»  _ . — 

ZMS  IliT(10ZKE,3.  .TRUTtGOSt 

ZM6 I'traOZMfcta.  .TRUE.GO&I 

ZM7  i:ir(iOZ.47,3..TRUE,607) 

- - ZM6  INra02M8,3.,TROE.G08l 

ZM9  INtdOZHS.O.  .TRUE.GOg) 

ZM13  INT(10ZM10.0..TRUEtGOlOI  

ZMIA  ItlT(10ZMlA,3.» 

CALCULATE. ACCEUCRAIION  OF  TIRE-tREAO 

20ZT1  (-FVU-  = 5U/sr-38o. 

20ZT2  (-FVT2-FS2l/>4T-386. 

23ZI3  <-Fvr3-FS3l/RT-386. 

20ZTI.  t-FVTi.-FS‘.l/,3T-386. 

23ZT5  (-FVT5-FS3I7HT-386. 

2aZT6_..  (-FVI6-FSSI/HI-386. 

20ZI7  (-FVI7-FS7)/9T-366. 

Z0ZT8  (-FWT9-FS31/HI-386.  

20ZI9  (-FVT9-‘^S9I/KT*386. 

20ZTi:  (-FVT1  J-F3l0t/M’^-385.  - 

TREAD  VJLOCIIIEj  AFTER  IMPACT 

lOZU- — LSMIG0,inZTL3.I).  T 

10ZT2  LSV((GQ2,l3ZT2O,0.l 

10ZT3  L3'H(GO3,13ZT3B,0.l 

10ZT4  LSMIGOA.IOZTAB.Q.I 

10ZT5  L3M(G05.iaZTS8,0.l 

IQZTb  LSH(GOo,nzT6B,IJ.) 

10ZI7 LSM(G07,l0ZT7B,a.J 

1DZT8  LSM(GD3,IOZr3B.O.I 

10ZT9  . LSH(G09,i0ZT93,#.» 

lOZTlO  LSM(G0lj,l0ZTA8,D.I 

- - INTEGRATE  ACCELERATIONS  AFTER  IMPACT  — - 

lOZTlB  INT(20ZT1,Z1, TRUE, GDI 

A.0ZT28-.-  INT<20Zr2,Z2,IRU£,G02l 

10ZT33  INH20ZT3,Z3,TRUE,G03I 

IOZThB  INT(2CZm,Zi»,TRUE,G04l  . 

1DZT5B  IIII(20ZT5.Z5,IRUE,GOE» 

1QZT6B  INT(2CZT&,Z6,TP.U£,G06) 

inZT7S  INI12D2T7,Z7,TRU£,G07I 

lozTea.  . iiinzozTs.za, true, GOBI 

10ZT98  INI (zozrg.zg, TRUE, C09I 

IQZTAB  I>IT(2DZI13,ZlC,Tr.UE,G011I 

TRACK  AND  STORE, STORES  AT  IMPACT 

zi  tas(vti,i:u,3.i  

Z2  TASIVT2,IN2, ).» 

--Z3 IAS(VT3,1N1,0.I 

74  TAS(VI4,IN4,3.) 

Z5  TA3(VT  5,I*»5,3.I  

Z6  TA3(9T4,IN6,C.» 

Z7  TAS(VT7,IN7,3,I  

Z8  TAS(VT9,IN3,3.I 

Z9  TAS(VT9.IN9.3.> 


ZIO 


TAS(VU),INiatO<l 


IMITIAt  IMOACI  velocity  OF  SPOKE 


VTl 
VTZ.  . 
VT3 
-VIi.— 
VT5 
VT6 
VT7 
VIA 
VT9 
-VTta- 
INl 
IN2 
IM3 
INL 
IN5 


IH6 

IN7 

.IN« 

IN9 

IMia 


nZKlA-9H»tOTHl»COS<THll 
iaZW14.9K*tOTHl»COS«TH2»- 
10ZKl*-?«»l0THl»COSITH3l 
-10ZmA-3H»u>IHl»C0SiTMi,» 
10ZKl»->H*l0THl»C0StTH5l 
10ZM1A-RM»10IH1»C0SITH6I 
10ZWl»-?W«lOfHt*COS<TH7» 
10ZHlA-KH«l.QrHt<C0S(THftl_ 
10ZW1A-YH»10THI»CQSCTH9| 
-lDZmA-RM»iOTHl»CQS(IHm 
LSMtGO, FALSE tTRUEl 

LSmCOZ, false, TRl»E> 

L$M(G03, false, TRUEI 
LSH(CQ,,FALS£.TRUEJ.... 
LSH(6C3. false, TRUE! 
LSMtGOS, FALSE, TRUCl 


LSH(G07, false, »RUE> 

LS,MG08, FALSE, TRUE3 

t.SH<G03, FALSE, TRUE) 
LSH(GCia, FALSE, TRUEI  . _ 
TREAO  OISPLACEMENT  AFTER  IMPACT 


_1U 
ZTZ 
ZT3 
ZTL 
, ZT5  _ 
ZTo 
ZI7._ 
2T6 
, 219 
zua 


INTU02Tl,3.,rRue,GOI, 
IMT  uczrz.a.  ,TRUE,GOZI 
IMmOZ^3,0.,TRUE,GO3> 
INTI10ZI4,3. ,TRUE,G04J 


INKlOZTS.a.  ,TRU£,G05) 

INTdOZTE.a.  ,TROE,G06l 

iKi(iozrr,a.  .iru£,gq7», 
INTdCZTS,  1.  ,TRUE,GO«l 

INTdDZTR.a.  ,TRUE,G09I 

IMT  dczm,c  . ,TRUH,GO10) 

NORMAL  FORCES  ON  IRE  AOS  

N1  -FS1*SIV(TH1* 

H2 -FSE»SINUH2Li 


-FS3*SIN<IH31 

-FSL^SIUdHM. 

-F3E*SIH»TH5) 

-FS6»SIN(IH6I_ 

-FS7*SIMTh7> 

„-FS«»SII(TH8Z 


-FS9»SIN«TH9I 
-FS10»SIN(TH10I 
CALCULATE  ACCEL! RATION  OF  SOIL  HASS 

1FS1-C3S»10ZS1»/MS 

(FS2-CSS»lOZSZ»/HS 
„<FS3-C53»10ZS3»/NS 
IFS4-CSS»1QZS4I/HS 
«FS5-CSr,»10ZS5>/HS 
lFSfc-CSS»lOZS6l/MS 
(FS7-CSS*10ZS7>/HS_  „ 

«FS8-CSS*lOZS«I/NS 
. . (F39-C3S»10ZS9I7MS 


Z02SI 

20ZS2 

20ZS3 

20ZS4 

20ZS5 

20256 

20257 
20ZS6 

. 202S9 
23ZS10 


(FS10-:SS’10ZSt«l/MS 
SOIL  FORCE  MUST  BE  COMPRESSIVE 


FSl 

FS2 

FS3 

._.FS4.. 

FSS 

FSo 

FS7 

FSa 

FS9 

FSlO 


FSHIFSlA.FSIAiQ. ,0.1 
FS,MFS2A,FS2A,0.,0,» 
FSH(FSIA,FS3A,0.,0.) 
-_FS«<FS4A,FS4A,0,.,OL,.L 


FSM<FS54,FSSA,0.,0.I 

FSNIFSSA,FS6A,0.,a.i 

FSmFS7A,FS7A,0.,fl,» 

FSNIFS9A,FSSA,0.,Q.t 

FSH(FS9A,FS9A,0.,0.» 

FSM<FstaA.FStaA.o..n.i 


ijKw;-'  H'<  ■-;jsw®«''r* 


2^10  tNmnzsto.a..T9UE.Goioi 

CALCULArt  TIR£  RQTAflOSAU  ACCtURAIION 

20TH1  <-TOR-i1l/IM 

CALCULATE  MOKEMS  DUE  TO  SPOKE  FRICTION  . 
TOR  T l»R*T ’•P2*T  T»R3*IA»RA*T5»RS*T0RI 

lORt T6*P.o«r2»S7*IS?-R#H9»Jl9*X10tR10 

CALCULATE  RAOTI 


..  ..  R RMt(ZH1-ZT1)/S1N(TH1J 

RZ  RN*(ZM2-2T2I/SIM«TH2J 

R3  RM»(ZHI-ZT3t/SIN(TH3) 

R4  PMf <ZMA-ZT4I /SlNTTHAt 

RS RW*(ZM5-ZJ51/SINTTHSI 

RS  RWA(ZV<3-ZT6I /SIN(TH6I 

...  RZ  . _ RHAIZHZ-ZT7T/SINITH7I 

RS  RR*(ZM3-ZT8(/SINTTHSI 

. - R9  RM*(ZH9-ZT9I/SIN(TH9I 

RIO  R«*<ZHlJ-Zri3)/SIN(TH10» 

lOTHl.  . INT.<20TH1,0NI 

ALL  SPOKES  START  9ASE0  OM  THl 

.THl  INTtlOTHl,a.ai!6fTRU£,COA» 

THl  RESETS  AT  MAX  ANGLE 

_ GOA  FSHUHl-Z.lZo. TRUE, TRUE. FALSE! 

ANGUS  CAN  NOT  3E  NcSATiVE  OR  ZERO 

._IH2 FSH(rH21,3,6J6,Q.6C6,TH211 

TH3  FSM(TH31,a.606,O.6C6,TH31> 

. .THL  FSH(THi»l, 0.636,0. 606, TH41I 

TH5  FSH(TH5l, 3.606,0, 606, TH51I 

_TH6  — FSH<TH61,3.616,0.b06,TH6ll 

TM7  FSN(TH71,3.6J6,0.606,TH71) 

„ THO FSN(THSl,3.6J6,a.606,TM8U 

TH9  FSKUHn, 1.636,0. 606, TH9II 

TH13  FSMTTIU. 1,3. 606,0, 606, THIOII  


TH3 
. -THL 
TH5 
_TH6 
TM7 

..  THO 

TH9 

THIO 


GET  LOGICAL  CONTROL  AWAY  FROM  ORIGIN 
GOZOi  FSM«(T-(l. <.1m)7<0NII, TRUE, TRUE, C020)  - -.  ..  .. 

G0381  FSN<  (T-(1.!.Ih1/C0H)), TRUE, TRUE, C0391 

GOnBI FSHT  (T-(l.<,14)/(OMM  ,T.xUE,TRUE,GO<taJ 

G0501  FSMI  (T-U.I,li,l/IOM»  I .TRUE, TRUE, GOER) 

G0691  FSNUT-(l.i,l<.>/<ON»), TRUE, TRUE,  G0691 

00761  FSH(  (T-(l.(,l<i>/(OHM  .TRUE, TRUE, C0731 

G08B1  F5H(  (T-ll. 41  i(»7(0HU  .TRUE. TRUE, G089» 

G0961  FSHT IT-(1.4l4»/(OM» ) , TRUE, TRUE, G093) 

jGOlOei  ...F5NT(T-(1. 4141/ (OMJ), TRUE, TRUE.GOUflJ 

ESTABLISH  START  FOR  OTHER  ANGLES 

GOZe  FSMITHl-C. 763, false, FALSE. TRUE! 

G038  FSM(THl-0. 923. false, FALSE, TRUE) 

G046  FSH(TH1-1. 177, FALSE, FALSE, TRUE) 

G05B  F3H(TH 1-1. 234, FALSE, false, TRUE! 

G06B  FSNTTHl-i. 391, FALSE, FALSEJRUEi 

G07S  F5H(THI-1. 54 3, false, false, TRUE) 

G088  FSH<TH1-1. 70S, false, FALSE. TRUEI  

G098  FSN (THl-1. 862, FALSE, FALSE. TRUEI 

G0103  FSHITHl-2. 01 3, FALSE, FALSE, TRUE)  

EQUATIONS  FOR  OTHER  ANGLES 

TH21 L5H(C0231,TH1-0.157,TH1»2,019) 

TH31  LSN(GOini,THl-0.3i4,THi41.862l 

TH41  LSM(G0481,THl-S.471,THlti.705) 

TH51  LSNIGOOll.THl'O.&ZS.THlTl.SLBI 

TH61  LSm006iU,THl-0,785,THl*1.391)  - . 

TH71  LSM(G07B1,TH1-0.942,TH1*1,234) 

I.H81 LSN(G03‘31,THl-1.099,THl*l.aZ7J 

TH91  LSH (GO  381, TH 1-1.266, THl AO. 920 1 

THIOI  L3H<C01:B1,TM1-1,4*3,THIAO,763)  . . _.. 

calculate  accfleraiion  of  wheel  hub  horizontally 

20XH  (OP»IOT-FORCEI/<HMaNP)  _ 

SUN  tangential  CONTRIBUTIONS 

TOT  T1»S1N(TH1)aT2»SIN(TH2IaT3»SIN(TH3MT4*STN(TH4UTOT1 


TOTl  T5*SIH<TM5(*r6»SINIIHb>*I7»SIN(TH?)*rS*SI><(TH8)*T0I2 

T0I2  19»SIN(IMll»tia»SINIIMlal 

92  SUH  SPOK<^  AXIAL  C09T9I3UTI0N 

FORCi  FXIl»C3J«THl»*fVT2»CCS«TH2> ♦FVT3»C0S«TH3) *F0RC1  ■ 

FORCl  FVT<.»CS3(TMAI»FVTS»CCSIIH5l*FWT6»C3S«THo»*FORC2 

F0RC2 FWI7»C0S(IH74»F«*»CCS(I«ai  ♦F-«J9»COS  (IK9I-*F0RC3 

F0PC3  FVTIO'COSITHIO) 

lOXM  Inn20XM.(/0»- 

XM  iNTdOXH.).* 

. ORAC  ((iVJ«RP»*20XR-0P--  _ — 

AORAG  (INT  lORAG,:.  I )/(T«'.0Ct> 

_l» CALCULAt£-IOTAL-SLIP 

TOSLP  «RH*10tHl-10XH»/<10THl*RR» 

OTHAX  . OTNIK  - 

OTHIN  .OOCl 

fimjt.Txs* 

PLoa,m,2r.3.zT5,zT7..ri9i 

PtO(T, AORAG, TOSLPI 


afc?S»ks,s£ 
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